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Abstract 

It  is  shown  how  to  convert  the  integral  equation  for  the 
pressure  distribution  on  a  circular  disc  in  an  axi symmetric  sound 
wave  into  a  singular  integral  equation  of  the  first  kind.   This 
singular  integral  equation  has  a  simple  kernel  and  can  be  trans- 
formed easily  into  an  integral  equation  of  the  second  kind  which 
is  especially  useful  at  high  frequencies.   Although  a  direct  iteration 
on  this  last  integral  equation  fails  an  indirect  method  is  devised 
which  ensures  that  each   iterate  is  of  lower  order  than  the  preceding 
at  high  frequencies.   The  form  of  the  general  term  in  the   iteration 
is  given  together  with  the  first  terms  of  its  asymptotic  behaviour. 
It  is  relatively  simple  to  estimate  the  error  caused  by  stopping 
at  any  particular  iterate  . 

Detailed  calculations  are  made  for  a  plane  wave  striking 
the  disc  at  normal  incidence. 

Introduction 

The  problem  of  the  diffraction  by  a  circular  disc  has 
been  investigated  intensely  during  the  past  decade,  primarily 
because  it  is  the  simplest  problem  in  which  the  diffracting  edge 
is  of  finite  length  and  therefore  serves  as  a  convenient  test  of 
methods  which  may  be  applicable  to  more  general  shapes.   In 
principle,  any  numerical  information  required  can  be  obtained 
from  the  solution  in  terms  of  spheroidal  functions  and  indeed  a 
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considerable  amount  of  such  information  was  derived  by  Bouwkamp  pf] 
many  years  ago.   In  recent  years  the  emphasis  has  been  on  attempts 
to  provide  analytical  information  which  can  be  used  to  supply 
backing  for  approximate  theories  of  more  general  shapes.   Apart 
from  one  method  employed  by  Chako  [2]  and  by  Nomura  and  Katsura 
[3]  in  which  an  aperture  distribution  in  terms  of  hypergeometric 
polynomials  is  assumed  the  main  approach  has  been  via  integral 
equations.   The  author  (Jones  [V] )  discovered  a  relatively 
simple  integral  equation  which  is  especially  suitable  for  finding 
results  quickly  at  low  frequencies.  Other  derivations  of  this 
integral  equation  have  been  given  by  W.I.  and  A.N.  Akhiezer  [p] , 
MacCamy   and  Heins  [6]  ,  Bazer  and  Brown  [f]  ,    Noble  [8]  ,  Collins 
[9]  ,    Bazer  and  Hochstadt  [lOj  ,  Williams  [llj  .   A  different  integral 
equation,  also  suitable  for  low  frequencies  though  less  simple  than 
the  preceding  one,  was  found  by  Magnus  [l2j  ;    it  has  been  derived 
by  other  methods  by  Noble  (loc.  cit.  ),    Lebedev  and  Skal'skaya  [13J 
and  has  recently  been  generalised  to  some  other  obstacles  by  Collins 

[iq . 

At  high  frequencies  the  standard  integral  equation  has 
been  solved  approximately  by  modifying  the  kernel  so  that  it  is 
amenable  to  the  Wiener-Hopf  technique  (Levine  fl51  ,    Levine  and  Wu 

[16]  ).   By  this  method  Levine  and  Wu  obtained  the  scattering  coefficient 

-  5/2 
as  far  as  0(o:      ),   where  a   is  the  product  of  the  wave  number  and 

radius  of  the  disc.   Their  method  was  used  subsequently  by  Seshradi 

and  Wu  [if]  for  the  diffraction  of  electromagnetic  waves. 
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Approximate  methods  based  on  physical  reasoning  have  been 
suggested  by  several  authors.   One  method  (Clemmow  [l8] ,  Millar  [19, 
20] )  uses  assumed  currents  on  the  screen.   A  theory  which  employs 
rays  as  the  basic  concept  has  been  given  by  Keller  [21,  22j  ,  Keller, 
Lewis   and  Seckler  [23],  Karp  and  Keller  [2^. 

Recently  a  very  general  method  of  tackling  the  asymptotic 
solution  of  diffraction  problems  has  been  proposed  by  Keller  and 
Buchal  [25] .   Their  method  is  applicable,  in  principle,  to  discs 
of  arbitrary  shape  and  it  is  therefore  important  to  confirm  that  it 
gives  the  correct  answer  for  the  circular  disc.   Certainly,  it  supplies 
a  scattering  coefficient  which  is  in  agreement  with  the  calculations  of 
Levine  and  Wu  as  far  as  they  go.   For  confirmation  of  higher  order  terms 
it  would  be  necessary  to  extend  the  analysis  of  Levine  and  Wu. 
However,  this  extension  would  not  be  simple  since  their  method  involves 
two  approximations  -  the  first  replaces  the  kernel  by  a  more  convenient 
one  and  the  second  consists  in  converting  the  resulting  integral  equation 
into  a  series  of  Wiener-Hopf  integral  equations  by  an  iterative  pro- 
cedure. 

It  seemed,  therefore,  that  a  new  approach  might  well  be 
profitable.   The  object  of  the  present  analysis  is  to  provide  a  new 
integral  equation  which  is  simple  and  especially  suitable  for  handling 
high  frequencies.   The  general  term  in  the  asymptotic  development  can 
be  obtained  and,  in  principle,  the  field  can  be  evaluated  to  any  desired 
order,  the  main  labour  involved  being  the  asymptotic  estimation  of 
certain  integrals. 


In  section  1  is  derived  a  singular  integral  equation  of  the 
first  kind  for  the  sound-soft  disc.   This  singular  integral  equation 
is  converted  to  an  integral  equation  of  the  second  kind  in  section  2  and 
the  particular  form  that  it  takes  for  an  incident  plane  wave  is  determined 
in  section  3.   Section  k   is  concerned  with  an  iterative  procedure  for 
solving  the  integral  equation  of  the  second  kind;  this  procedure  is 
especially  appropriate  at  high  frequencies.   The  resulting  distribution 
on  the  disc  is  found  in  section  5  and  the  high  frequency  scattering 
coefficient  in  section  6. 

In  section  7  "the  corresponding  integral  equations  for  the 
sound- hard  disc  are  derived  and  section  8  is  devoted  to  the  appropriate 
iterative  procedure.   The  distribution  on  the  disc  and  scattering 
coefficient  are  determined  in  sections  9  and  10  respectively. 

It  is  shown  in  section  11  that  the  solutions  of  the 
integral  equations  of  the  second  kind  described  above  are  unique  and 
that  the  iterative  procedure  does  provide  an  asymptotic  development 
at  high  frequencies. 

The  derivation  of  certain  results  has  been  carried  out  in 
appendices  .  References  to  equations  in  appendices  are  written,  for 
example,  (Bl),  meaning  equation  (l)  of  Appendix  B. 
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1.  Derivation  of  the  Integral  Equation 


We  consider  the  problem  of  the  diffraction  of  a  small- 
amplitude  harmonic  sound  wave  by  a  sound-soft  circular  disc  of  radius 
a.   The  boundary  condition  is  that  the  field  vanishes  on  the  disc. 
The  case  of  the  sound- hard  or  rigid  disc  where  the  normal  derivative 
of  the  field  vanishes  on  the  disc  will  be  treated  in  later  sections. 

Let  r,9  and  z  be  cylindrical  polar  coordinates,  and 
normalize  all  lengths  by  division  by  a  so  that  the  disc  occupies 
0  <  r  <     1,  z  =  0.   Let  the  incident  field  be  independent  of  cp  and 


be  represented  by  u  (r,z)  at  (r,cp,z),  the  time  dependence 


ia>t 


being  understood  and  omitted  throughout.   The  total  field  u(R),  at 
the  point  whose  position  vector  with  respect  to  the  orgin  is  R,  is 
then  independent  of  cp  and  may  also  be  written  u(r,z). 
The  field  satisfies  the  equation 

(V  2  +  a2)u  =  0 

where  a  is  the  product  of  the  wave  number  k  and  the  radius  a. 

Throughout  a  will  be  taken  as  a  positive  real  number.   The  total 

field  at  the  point  R  can  always  be  written  in  the  form 

~du 


t(R)  =uQ(r,z)  -^ 


3z\ 


?//(R,R  )dS 


since  u  =  0  on  the  disc.   Here  S  is  the  unit  circle  on  z  =  0,  R, 


-1 


is  a  point   of  S,   ^(R,^)   =  e~lal-       -l'/    |R  -   R±\    and 


du 
3z\ 


du(r1,z1) 
"" 5z~ 


o\i(r1,z1) 


-5z— Jz1=  -  o 


Applying  the  boundary  condition  u  =  0  we   obtain  as  an  integral  equation 
to  determine  f(r    )    (   =    [du/dz  ]  ) 

(1)  uo(r,0)   -  \^J    f(rxv  ^(R7R1)dS  =   0 

where  now  R  is  a  point  of  S.   As  a  consequence 

2    2  -i  l/2 

exp[}-ia  {r  +  r,  -  2rr-[_  cos(cp-cp1))  ' 

(2)  ^(r,r  )  =  — -g — 2 ; — w~m 

(r    +  r     -    2rr     cosO-Cf^)}      ' 

When  7p  is   given  by   (2)   it   has  been  shown    ( Jones    [k\  )  that 

min(r,rn  )  ,    2     2 .1/2  ,      2        2 ,l/2 

/2rt  p        '   1  cos  <x(.r  -x    )  '      cos  a(.r,    -   x    J 

^i  =  y         ,  2  2,1/2 ,  2  2.1/2 ^ 
uo                       (r  -x    J  '       (r,    -   x    j 

(3) 


r  ,   2        2,1/2    .         /      2        2,1/; 

p    cos   al.r   -   x    )        sm  a(.r^    -   x    J 

"Hi  ,,2      ^l/2    ^        _    ^1/2 


(r  -  x  )  '       (r,      -   x   ) 


For  convenience  a  derivation  is   given  in  Appendix  A. 

When    (3)   is    substituted  in    (l)  the  result   is 

r  /   2        2x1/2    (    1    _/      x  /      2        2,1/2 

n      cos  a(r  -   x    )   '       I     P    f^Jr-jCOs   a^  -   x    ) 

(r'0)-'"/     ~T2 2072—^,/    —^ — 2TJ2 dri 

b  (r     -  x   )  '  I  "x  (r,    -   x    )  ' 


n 


1   W       \  •         f      2        2^/2 

/f  (,r    Jr     sin  a(r     -  x    J 
.      2        2,1/2 dr      V  to    =  0 
o               (^  -   x    ) 

2  2  2 

for  r  <  1.      The  change  of  variables  r     =  y,    x     =  z,   r     =  t  converts  this  to 

,1/2,1      f  ^j±^2  {[fit1'2)  cos  a(t-z)1/2   ^ 

^     f^sin^t-z)1/2      Jdz  =0 

(t  -   z)1^  zV2 
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for  y  <  1.   The  solution  of  an  integral  equation  of  the  form 

(5)  fl(y)  -J    hx(.)    CMa(^2       «. 

^o  Cy-z; 

is  known    (Jones    Q+]  )   to  he 

i      (*)        1     i-    f  f    fv)      cosh  a(z-y)  , 

Hence  the   solution  of   (h)   is 
(6) 

/*     f(tl/2)   cos  a(t-z)1/2  dt   .    At1/2)   Bin  ,(t-.^2  dt 
^  (t-z)1/2  Jo  (t-z)1/2 


z  l/? 

1/2     d         f         /    1/2      ...     cosh  a(z-y)   7  , 

s      di   J   uo(y    '  0)  — ; — ^72—    dy 


(z-y)J 


valid  for  0  <  z    <  1. 

The  integral  equation  previously  described  in  the  intro- 
duction as  being  especially  useful  at  low  frequencies  is  derived 
from  (6)  by  applying  the  inverse  operator  of  the  first  term  of  the 
left  hand  side.   The  inverse  operator  is  easily  found  by  changing 
t  and  z  to  1  -  t  and  1  -  z  so  as  to  obtain  an  integral  equation  of 
the  type  in  (5).   To  find  an  integral  equation  suitable  for  high 
frequencies  (large  a)  rewrite  (6)  as 


(7,r  f^-^t-o^    rzf(tl/2)  cos  a(t.,)f  at 

•'  (t-z)1/2  Jo  (t-z)1'2 


z  l/2 

1    V2     d      /        f   1/2     n)      :osh  g(z-y)  7 

= 4z       dTj  uo(y     »  0)   — ; — a72~~    dy' 


1/2  1/2 

For  definiteness  (t-z)    means  i(z-t)    when  z  >  t.   It  is 
proved  in  Appendix  B  that,  for  the  x  and  t  involved, 


/£ 


/  ,    a/2-i   -ia(t-z)1'2 
os   ■ja(x-z)  '     )    e 

.1/2   ,  '   a/2 


(x-z)     (t-z] 

d/2   1/2 

t  '  +  x  '     -lay 

-     dy  -  itiH(x-t ) 

1/2   1/2  y 
t  '  -  x 


where  H(x)  is  the  Heaviside  unit  function  and  the  principal  value  of 
the  integral  is  to  be  taken.   In  fact,  we  shall  always  use  the 
convention  that  a  principal  value  is  implied  when  the  integrand 
contains  a  simple  pole  on  the  interval  of  integration.   Another 
relation  proved  in  Appendix  B  is 

(9)  f      C°S  {a(x-zffi  °°b  (a(t-z)1  2)   dz  =  _rti  sgn(x_t). 
Jt         (x-z)1/2(t-z)1/2 

Multiply  both  sides  of  (7)  by  cos  (a(x-z)  '  }   /(x-z) 
and  integrate  with  respect  to  z  from  0  to  x.   On  account  of  (8) 
and  (9)  equation  (7)  becomes 
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(io)     !  r    t1/2^1/2 


f    fCt1/2)/    f  e"iaydy-   «lH(x-t)) 

^  \Jl/2         .1/2     y 


dt 


t1'2   -  x1/2 


.1/2 


+  /  f(t  '     )   rtl  sgn(x-t)  dt 
o 


&  (x-z)  '  o  (z-y) 


Since  the  left  hand  side  of  (10)  is 


[fit1'2)    f  ,      f^dydt 

J  Jl/2        1/2     y 


t^2-    XV2 


a  derivative  of  (10)  with  respect  to  x  supplies 

*  ^J^A^)      ^  *1/2) 


at 


dx  J       (x-z)1/2      to-i   °  (z-y)1'" 

o 

where,  as  already  stated,  a  principal  value  is  implied  if  necessary. 

The 
(11) 


2  2 

The  change  of  variables  x   =  v     and  t    =  w     now  leads  to 


>                   -ia(w+v)  .    ,        »  /    \ 

vf(w)   (2 +  e-la(w'v)   }    dw   =  F(o)(v) 

W    +  V  w-v 


e  2  z 

f(o)(T)  .  2,f  ^        f  ta,  .^  a_^o(yl/2)0)  ^%£2dy  a, 

o  (v  -zj  '  t)  \z-y j 


where  ..2 

(12) 
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The  integral  equation  (ll)  holds  for  0  <  v  <  1. 

It  is  easy  to  see  from  the  structure  of  (ll)  why  it  may 
be  expected  to  he  especially  helpful  at  high  frequencies.   For 
large  a  the  dominant  contribution  will  come  from  the  neighbourhood  of 
the  principal  value ,  corresponding  to  the  physical  situation  that 
to  a  first  approximation  the  distribution  on  the  disc  is  determined 
locally  by  the  incident  field.   Furthermore  the  most  significant 
values  will  tend  to  occur  near  an  edge  v~l  since  there  the 
neighbourhood  of  the  principal  value  is  about  to  move  out  of  the 
interval  of  integration.   (The  w  in  the  numerator  prevents  v  =  0 
being  similarly  important).   The  first  term  of  the  integrand 
represents  the  effect  of  secondary  radiation  and  shows,  in  particular, 
that  for a  point  near  the  edge  the  main  secondary  waves  come  from 
the  diametrically  opposite  point.   However,  it  should  be  noticed 
that  as  v  approaches  0  both  terms  in  (ll)  are,  roughly  speaking,  of 
the  same  order  of  significance.   In  essence  this  corresponds  to  the 
fact  that  every  point  of  the  edge  is  equally  important  with  regard 
to  producing  a  field  at  the  centre  of  the  disc. 

In  the  low  frequency  case  when  a  is  small  both  terms  on 
the  left-hand  side  of(ll)  are  of  equal  importance.   Some  consideration 
of  the  low  frequency  problem  will  be  found  in  Appendix  C. 

2.   Conversion  to  an  Integral  Equation  of  the  Second  Kind 
There  are  several  ways  in  which  (ll)  can  be  transformed  into 
an  integral  equation  of  the  second  kind.   In  view  of  the  remarks 
made  at  the  end  of  the  preceding  section  concerning  the  relative 
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significance of  the  two  parts  of  the  kernel  at  high  frequencies  we 

put        1 

!,(w-v) 


P    /  »  -ialw-vj 

<13>     J    ^       *  "  *>■ 

o 


Expressed  in  the  form 


1  __.o/__\_-iaw 

dw  =  F(v)e 


wf (w)e  "    j„   T7i/'.^_-i<xv 


r 

/     w-v 
Jo 


(13)  is  a  singular  integral  equation  of  well-known  type  (see,  for  example, 
Mikhlin  [26] )  with  a  general  solution 


r  1/2  r-,       \l/2 

„i    \   -law  1         /  v   (,1-vJ    w  -,  -iav       C 

wf(w)e     =  -   2  l/2     l/2  J  —7^ F(v)e     dv  +    j                 j 

Jt  W  '   (1-Wj  '    T>  W  '   (1-WJ  ' 


where  C  is  a  constant.   C  is  determined  by  applying  a  condition  at 

either  w  =  0  or  w  =  1.   Since  no  edge  condition  is  necessary  to 

ensure  uniqueness  (Jones  [27] )  for  the  sound-soft  disc  we  impose 

the  condition  that  f(w)^and  therefore  wf(w),  is  bounded  at  w  =  0. 

This  leads  to 

iaw        p   ,_,  _\l/2  „/__\    -iav 


(Ik)  ffv)         e^W        f   (l-v)1^   F(v) 

(U)      f(w)  =  "    2  1/2, n   ,1/2  J   ~T72     v~T 

7t  W  '   (1-W  )  '      O    V  ' 


S     dv. 

(i-wT 


By  substituting  for  f(w)  from  (ik)   we  find 

1  ,  /„  1 


-ia(w+v)  p      -iav         l/2  p  /.    ,  \l/2      ,    .   -iat 

~/    \  e  ,  1      /     e  w  '  /  ll-t;  F(t  je  ,.     , 

wf(wj  — —  dw  = _   /      — - —  , /o   /  — T75 ~— —         dt   dw 

w+v  2  J        w+v  /,       w2J  .±2.  t   -  w 

rt     o  (1-wj  o         t 
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11/2  1  ,  /o 

■?    4  T1^  F(t)e    4  ^/2  ^^ 


on   inverting  the   order   of  integration.      The  inner  integral  may  he 

easily  evaluated  hy  introducing  the  variable  x  defined  by 

2  2 

w  =  x  /(l+x    ).      We  thus   find,   remembering  that  a  principal  value 

is   involved, 


(An  alternative  derivation  will  be   found  in  Appendix  Ej    see    (E  Ih  ) ) 
Hence 


/-ia(w+  v)                 -iav     l/2 
wf(w)      ■ dw   =  -  T-r 

w+  i  n  (±+      a/2 


(l+v)J 
(16) 

X    f    d-t)l/2   E(t)e-iat      dt. 

t   '     (t+v) 


Since,    from   (ll)  and   (13), 


(17)  F(v)    =  F(o)(v)   -J    wf(w) 


w+v 
ro 


-ia(w+v) 

dw 


we  obtain  as  an  integral  equation  for  F 


(!8)F(v)    =  F(o)(v)   -   ^"!t     l(l-)t/2        F(t)e"iat   dt      (0  <  v  <  l). 


n(i+vr/c  ^o  tx/^(t+v) 
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When  (l8)  has  been  solved  for  F  the  function  f,  and  hence  the 
whole  field,  is  known  from  (lh). 

Equation  (l8)  differs  from  (ll)  in  being  an  integral 
equation  of  the  second  kind.   Its  kernel  is  perhaps  not  quite 
so  simple  looking  but  it  is  no  longer  singular  and,  as  we  shall 
see  later,  the  appearance  of  the  kernel  can  be  somewhat  simplified 

1/0        ~]   / o 

by  taking  (l+v)   F(v)/v    e     as  a  new  unknown. 

It  is  also  possible  to  form  an  integral  equation  of  the 
second  kind  for  f  by  substituting  for  F  from  (17)  into  (lU). 
The  result  is 

1 
law  p      /   -iav   /  ., 

f M  -  -  e21/£,     )1/2  J  <^>1/2  hr-  tfo)v 

Jt  w  '  (1-w)  '        uo 

-ia(t+v) 

tf(t)  2-——      dt)  dv. 
t+v 

However  it  does  not  seem  possible  to  simplify  the  kernel,  even  to 
the  extent  of  eliminating  one  of  the  integrations.   Consequently, 
the  kernel  is  considerably  more  complicated  than  that  of  (l8)  and 
approximation  by  iteration  would  involve  a  double  integration  for 
each  iterate.   In  fact,  the  form  of  the  first  term  for  an  incident 
plane  wave  is  such  that  there  does  not  seem  to  be  a  uniformly  valid 
representation  any  simpler  than  the  original.   It  is  therefore 
advantageous  to  operate  with  (l8). 
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Both  (ll)  and  (l8)  seem  to  be  new.   The  integral 
equation  derived  by  Magnus  [l2j  bears  some  resemblance  to  (l8) 
but  is  in  fact  more  complicated. 

3.   The  Incident  Plane  Wavt 

In  order  to  see  how  a  detailed  solution  of  (l8)  goes 
consider  the  special  case  of  an  incident  plane  wave  in  which 
u  (r,z)  =  e    .   Then,  from  (12), 


F(o)(v)  =  2 


d 

n  t    (   2    .1/2-,    ,    1/2 

d    /  :os{q(,Y  -  z  J    (■cosh  qz 

av  1  r  2   ,1/2 


dz 


(v  -  z)" 


=  2 


dv 


2   •   /  I    2    \l/2\ 
-  -  sin  i  a(.v  -   z)        (    cosh 


as 


1/2 


-1  o 


2^     f    sin{a(v2-z)1/2} 


sinhaz 


1/2 


'W 


dz 


after  an  integration  by  parts.      Hence 

2 


F^°'(v)    =  h  cos  av   +  2av 


I 


r  /  2     aM     .  ,      1/2 

s[a(v-    z  J   '     j-      sinh    az 

,  2     a/2     1/2 
(v   -z  J  z 


dz. 


-1^- 


Since  (sinh  az  '  )/z  '   can  be  regarded  as<sina(-z)    >/(-z) 
(B3),  with  x  =  v  and  t  =  0,  shows  that 


y 

at 


dt 


F   (v)  =  k   cos  av  +  2av  /   — - 
=  h   cos  av  +  4av  Si(av) 
where  Si(x)  =  /   — - —  dt.   An  alternative  way  of  writing  Si  will 

/  COS    "t 

be  found  valuable   in  the   following.      Since    /        — - —  dt   =  0, 

■it  ,  p        -i(t+x) 


,       p      -it  ,  p        -K.t+xJ 

(20)  SiUJ.ii/    f       at.ii/-i_—     dt. 

-x  -2x 

Consequently  o  .    / ,  .    \ 

(21)  Si(av)    =  |  ij 


t+v 

Now  define  a   new  unknown   F   (v)  by  means   of 

00     -ia(t+v) 

(22)  F(v)    =  2eiaV   +  2iav   /        ^— - dt    +  F   (v). 

J-2v 

The   equation   satisfied  by  F   (v)   is,    from   (l8),    (19 )  and    (21). 

(23)  00  ,         .  , 

n      -ia(t+v)  l/2  -iav 

iccv       ,-,_.    ..    /      e  „     v  '    e  _        x 


F(v)   =  2e-iav    -   2iav  ?— dt-   * 2_ 

1  'I  t+V  n(l+v)1/2 


1         w^    .       r 

=  0 

o       w  '     (w+v)  |  -2w 


for  0  <  v  <  1. 
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Now 
(24) 


J    ^72; — :  dw  =  *  <    1/2     - x 

*"b  w  '     (w+v)  1     v  ' 


for  the  range   of  v  involved,    as   can  be   seen  by  the    substitution 

w   =  x  /(l+x    )   for  example    (or   from   (E  10)  as  t  -**>). 

Also 

r/n       \l/2      1/2      -iocw    p        -ia(,t    +  wj 
(1-w)  '      w  '       e  /        e 


r>     /-,       sL  d     Id      -law    r>        -ia  t 
/      (1-w)  '      w  '       e  /        e 

J  (W       +     v)  Jp.r     "  + 

^^  1  /  /  * 

fd-w)^2    r   e" 


dt   dw 

2w 


t   -  w 
~'o 


dt   dw 


2  . .        2 
by  replacing  t  by  t-2w.      Now  the    substitution  w   =  x  /(l+x    )    shows  that 

f   (l  -  -)1/2  -1/2  aw  =  -S-  /^e(i^)V2 

Jo    (w+v)(w-t)  t4v 

(26) 

+  tl/2(t-l)l/2  H(t-l) 

for  v  >  0,  t  >  0  (Ell).   Hence  an  inversion  of  the  order  of 
integration  in  (25)  leads  to 

p    .        ,1/2  1/2     -iocw    p        -ia(t+w) 
(1-w)  'we  /       e 


p    ,        ,1/2  1/2     -iaw    n        -ia(,t 
/      (1-w)       w  '      e  /        e 

4>  w+v  ^-2w  t    + 

00 

in   (  -ia  .  y  l/2h-L   a/2    P    e~lat      _. 

=  —   (e        -1)   -   rtv  '     (1+v)   '       /      — - —     dt 


-at 


dt. 
t   +  v 
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Employing  this  result  and  (2U)  in  (23)  we  obtain 


_  1/2  -ia(v+l)     Q.   1/2 
„   /    \       2v  '  e            2iav  '      -lav   v 
F  (v)  =  ^75 r/?   e      X 

1         (l+v)1/2         (l+v)1/2 


iav 
v  '  e 


.1/2/.  .  ^1/2*^   .  ,        1/2 
:  /  (t-1)  '   V^-Kxt  A+        v  ' 

:(l+v) 


Bill<i-^-i=rS-    •~«-t^37.    " 


/ 


1(l-w)1/2F1(w) 


-law 
e 


w- 


(w+v) 


dw  . 


A  rather  simpler  equation  can  be  derived  by  deforming  the  contour 


of  the  first  integral  into  the  line  joining  1  to  1-1  °°.   The 

o-   +-.  >  - 
2 


l/2  l/2  1 

branches  of  t    and  (t-l)    are  specified  by  taking  0  >  arg  t  >  -—  Jt 


and  0  >  arg  (t-l)  >  -  -  it.   Then 

2aVl/2     -^i-ia(l+v) 

F  (v)  =  =-,   e  X 

1      (l+v)1/2 
(28) 

'   ,1/2  1/2     .  .       1/2  -iav  ^(l-w)l/2Fl(w)e-iaW 
(1-it)  '  t      -lat  ,.    v  '  e     / 1 -1 

1  +  v  -  it   e    dt  "  h+   ,1/2  J    TfT— dw- 

•»0  it(l+v)  '  vo  w    (w+v) 

with  the  understanding  that  0  >  arg(l-it)  >  -  -  it.      Putting 


2 


(29)      Ffv)  =-2-^/2   e"iaVl(v) 
(l+v)1/2 

we  have 


■  ia-^i  J    U-itJ  '  t e"atdt 


1  .  n,.     ..   »l/2.l/2 

i-^Jti  /  (1-it )  /  t  ' 

J        1  +  v  -  it 

o 

(30)      1    r    (1    ^2  w  \  -2i^ 

.  i  /    d-w)  ^   I(w)e dw> 

*    Jo        (l+w)1/2  (w+v) 

The  approximate  solution  of  this  when  a  »1  is  considered  in  the  next  section 
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One  might  be  tempted  to  solve  (l8)  by  iteration.   However, 
if  we  write  (l8)  as  F  =  F   -  KF  and  work  in  a  Banach  space  with 


F| |    max  |f(x)|  we  have 

2 


0<x<]  ± 


. ,  v  f  /i-ta/2    dt 

(Xv<L   rt(l+v)  '  •   o 


max 
CKv<L 


{-^1/£}- 


so  that  we  cannot  say  that  the  iteration  process  will  converge. 
Nevertheless  the  process  may  give  an  asymptotic  development  and 
this  seems  not  unlikely  in  view  of  the  equivalent  form  (30).   This 
possibility  will  be  examined  in  the  next  section. 


\.   1  approximate  Solution  for  Large  a 
Write  (30)  in  the  form 


,  ,      ,  ,   1   f  ,l-w,l/2   l(w)e"2iaW 


+v 
o 


where 


00  11 

(32)  hl(v)   =  2ae-ia-7Ti/       d-it)12t12  e-atdt_ 


T)  1+v   -    it 
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A   formal  expansion  by  iteration  immediately  comes  to  mind  and  we 
examine  the   scheme 


(33)  ^(v)   =  h1(v), 


o 


T        /    \   -2iaw 
I      n(wje 

n-1  , 
; dw. 

w+v 


The  presence  of  the  factor   e  suggests  that  the   iterated  terms 

will  be   of  diminishing  order  of  magnitude    for  large   a.      However,    a 
check  reveals  that  this   is   not   so.      To   see  this,    note  that  as  a  -*  °° 


(35)  h]_(v) 


1/2     -ia-4rti 
n  '      e 


~J^Vi 


+v) 


for  all  v   in    [0,l]  .      Prom  this  we   observe  that    I?,      and  I    ,    ... 
are  of  increasing  order  of  magnitude   near  v  -   0.      It   is  therefore 
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necessary  to  adopt  a  less  direct  method  of  iteration,  which 
eliminates  the  difficulty  at  the  origin.   Note  that  the  difficulty 
cannot  be  evaded  by  subtracting  a  constant  from  I;  this  merely 
alters  the  order  of  magnitude  for  each  iterate  but  they  still 
increase.   To  avoid  the  trouble  at  the  origin  put 


v  J        t-v  2 


(36)  Kv)  =  (^r^  /    =±=-  e  dt  +iQ(v) 


where   H   (t )  will  be   specified   shortly  and  I      is  a   new  unknown 
(and  not   given  by   (33)  and    (3^)  which  will  be   ignored  from  this 
point   onwards).      Since 

r^\  f/l-w\  <lw  n       frt-lC2u,+_    -,  n    ,1+vn1/2~\ 

(3T)       ^  ~     (^)(w-t)   =  —  <  (-r}    H^-^-(— >    j 


for  v  >  0  and  t  >  0  (the  derivation  is  along  the  same  lines  as 
that  for  (26)  or,  alternatively,  as  in  Appendix  E  (ELO))  we  have 

-2i<x(t-v) 


^w-yvj-tis^jv^f 


t-v 


-21CCA  P    +    1    1/0     H   (t)   -2iat 


.  1  /'    (k^)i 

Now  choose   IL  (t )    so  that 


x 
1  T    /    n   -2iccw 

,/2        Z2(v)e 


w   +  v 


dw. 


[     -2ia(t-v)  -2iat 


t-v  t    +  v     f  1 
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or,  equivalently,  so  that 


ia(t-v)    -ia(t+v) 

+^T-T 

-  v        t  +  v 


=  W     hl(v)e     • 

This  equation  is  of  the  same  form  as  (Dl)  and  is  satisfied,  according 
to  (D15),  hy  taking 

00 
(39)      ^(y)  =  Tt^ IQ   {Tk)lf\   (x)e_iC"  M(x'y)  dx 

where 

hxy   h/2)(cc  x)  J  (ay)  -  %2H  (2)(ax)J  (ay) 
(kO)       M(x,y)  =  i ^-^ 2 • 

x  -  y 

With  this  choice  of  H   (38)  reduces  to 

1     n  /0   _  /  \  -2iaw 

r  1   h  f  1   l  /V-M  7   T2(w)e 

I2(v)  =  h2(v)  -  -  j  (— )     -^j^- 


dw 


where  0.  4- 

r  .   .  1/2  H  (t )e-2lot 

In  view  of  the  resemblance  between  (32)  and  the  formula 
for  h (v)  we  now  formulate  the  iteration  scheme 

(41)   I  (v)  =  (i^)1/2  f    jfl   e-2ia(t-v)  ^    (n  >  1} 

00 
(1,2)   Hn(y)  =  g  ^^^  C^)1/2  hn(x)e-iaX  M(x,y)  dx,   (n  >  l) 
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?  u   (a-  \   -2iat 

P       4-  i  i  /o   H   -■  UJe 


and  K  is  given  by  (32).   We  must  now  verify  that  the  scheme  has 
the  desired  properties. 

Assume  that  h  (v)  is  regular  in  the  portion  R(v)  >  -  1 

n 

of  the   complex  v-plane   and  that    h  (v)   is  bounded  as    [v]    -»  °°  in 
that  region;   this   is   certainly  true   of    h(v).      Now  provide  the 
integral  in    (k2)  with  a   complete   contour   from  0  to  °°  by  adding 
an  indentation  below  the  pole   x   =  y;   to  keep  H     unaltered  the 
contribution  of  the  indentation  must  be   subtracted.      The  dominant 

behaviour  of  the  Hankel  functions  at   infinity  in  0  >  arg  x  >  —z  it 

1  / 
is   given  by  a   constant  multiple   of  e  /x  '     .      Therefore,    when 

account   is  taken  of  the  assumed  properties  of   h  ,   the   complete 

contour   just   constructed  can  be  deformed  into  the   negative  imaginary 

axis.      Hence  1   .  . . 

1/2  9     ^1+1CV  l/2  h   (-ix)e-0* 


(hh) 


PlV  nJJ  2  Jvl-ix'  2     2 

J  it  Xd  x   +y 

i  xyK   (ccx)Jo(ay)   -   y  K^ax^Cay)  i    dx 


where  K   (z)   is  the  modified  Bessel  function  related  to  the   Hankel 
v  l     .  i       . 

(?)        "2rtl  ?VTtl  l/2 

function  by  itiH^      ;(ze  )  =  -2e         K   (z).      The  value  of   (l-ix)  ' 

is   specified,    as   in  the  preceding  section,   by  0  >  arg(l-ix)  >  -  -  it. 

To  evaluate    (hk)  asymptotically  use  the  integral  representations 

■xt„  l/2 


00  00 

K°(x)  "I  ft?^  - Ki(x)  =  xiVxV- 1,1/_  at 
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Then  the  term  involving  K     in   (hk)  can  be  written 

2a2  7  Iti+iay  T-   o        i/o    P.   ^5    .l/o      hj-ix) 


V,  o  X      + 


X 


"2  y  e  VGy;J    ^   "lj        J    kl^x'  "2 2 

it  ~                            o                            x      +  y 

-ax(l+t) 

e                     dx  dt . 


The  inner  integral  is  a  Laplace  integral  provided  that  h  does 
not  depend  exponentially  on  a  (which  is  certainly  true  of  h  ), 
Thus,    as  a  -»  °°, 

n       5     1/2  h  (-ix)           ,       , 

/    /    :        \  '  n -ax(l+tj  h    /mv 

J    (IITx)  "T-2-e                      dx-hn(0)X 

o  x  +y 

(^5)     „ 


I? 


— (1+t) 


.   2 
o      x   +y 


dx 


since  t  >  1.      The   integral  on  the  right   of   (V?)   is  a  Gilbert   integral 
and  occurred  in  the  theory  of  diffraction  as   early  as  1863    [25]  ; 
it  is   expressible   in  terms   of  Lommel  functions.      However,    it  will 

mostly  concern  us   only  when  y  >  1.      When  y  >  1  the  asymptotic 

2  2  2 

behaviour   of    (V?)   is   easily  found  by  replacing  x      +  y     by  y   . 

Hence 

/V      ,1/2      V-1*'     -ox(l+th  15"1/2  h»(0)  +0(hn      , 


2 
Since  the   substitution  t   =  2sec      0-1  gives 
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1 
2  * 


/  C 

/        ^   ~1^>    ,  dt    =  2"1'2   /      sin2    9  cos    e  a<9   =  rt/2 


9/2 


(t+1) 


the  term  involving  K  in  (kk)   supplies  a  contribution 


l5hn(0)e  7ti+lay  JQ(ay) 
'     1/2  3/2   .13/2 


+0(-S/2)   JQ(ay). 
a 


In  fact   this   is  valid  not   only  when  y  >  1  but  when  ay  >  >  1. 

With  regard  to  the  term  containing  K   ,    it   can  be  written 


2a  _2 


1     . 

-rrti+iay 


2_    ,  1-1/2 


^y\Me  I    <t<-  lP<«/    <£;> 


1/2  h    (-ix) 


2        2 
x  +  y 


X 


-ax(l+t) 

e  dx  dt , 


By  a    similar  argument   to  that  which  has    just  been  employed 


1/2      h  (-ix) 


/u d  n  \-\ 

(—   )  ^~ 

^l-ix;  2 

i_>  X       + 

■I-j 


-ox(l+t) 
e  dx 


"o     y 


hn(0)   +  i*J|  hn(0)_  1^(0) 


-ax(l+t )  /  n      \ 

e  dx  +  0(— /gj 

a 


1 
2 


2y2a3/2(l+t)3/2 


31 1 


V°)   +     2afe7     t  hn(0)  "  hn(0) 


+  o(-£/2). 

a 


The   change  of  variable  t-1   =  z(t+l)    shows  that 


a 

/ 


dt 


(n-i):,1/2 


[  {t-^'Ht^r  :  (n_1);2n-| 
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Hence  the  contribution  of  the  K  term  is 

o 


J-^ay)        e 
«1/2al/225/2 


ni+iay 


h    (0)   +2i      K    (0)_h    (0)l 

n  16a    )2  n 


a 


Combining  the  two  results  we  obtain 


Hn(y) 


;<i%>1/2  "„w  +  77577^75 


■riti+iay 


9i      J  1 


h   (0)+£_  <ih   (0)-h   (0) 
n  16a    \  2  n  n 


(h6) 


15-h  (0)e 

n 


1      .\. 

■jj-  rti+iay 


JQ(ay) 


1/2  3/2     0l3/2 


h 
+  0(3/2 )   JQ(ay) 
a 


when  ay»l.      For  ay»l  the  Bessel  functions  can  be  replaced  by 
their  asymptotic   expansions,   namely 


JQ(ay)   =   (^)1/2  cos(ay-^)   +  0(a  2), 

Ji(ay)  =  (^)1/2{cos(ay-i  n)-  &v  sin(ay-!rt)) +  °(a 


2). 


■2k- 


Then,    for  ay  >  >  1, 


"nW  -  &V2  V*>  "  -a/2    <^»V„ 

Jtay  ' 


where 


We  can  now  see  from  (U3)  (with  n  +  1  for  n)  that  h 

n+1 

has  those  properties  of  h  that  were  used  in  deducing  (kj). 

Hence  there  is  a  corresponding  formula  for  H  , .   Since  h 

n+1  1 

has  these  properties  it  follows,  by  induction,  that  (h'J )   is 

valid  for  all  n.   Furthermore,  the  dominant  contribution  of  (^7) 

to  h   _,(v)  for  large  a  is  0(h   a  ).   Thus,  in  the  relevant 
n+1  n,o 

range,  h    (v)  is  smaller  than  h  (v)  by  a  factor  of  0(cc  ) 

and  hence,  by  (h-l),    H  ,n  is  smaller  than  H  by  a  similar  factor. 

Consequently,  the  iteration  scheme  is  appropriate  for  large  a. 

y.   Determination  of  the  Distribution  on  the  Disc 


Denote  the  contribution  of  F(v)-F  (v)  to  f  by  f  and 
3ntri 
(29), 


the  contribution  of  I  to  f  by  f  .   Then,  from  (lk) ,    (22)  and 
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„    iccw  n  ,         >l/2   - 

Ww1    _  2e /    (1-vJ        e 

fo(v}  -  "    2  1/2  .     a/2    J  -172; — : 


r(l-v)l/2e-laY 

) 

(48) 

\  J-ax^  .         j  -ia(t+v) 

i.  ^-2v        t    +  v 


dt  I    dv, 


fa-v//2 


iccw  p 

fn(w)    =  1  1/2,         ,1/2    Jk 
it  w  '     (l-w)  '       w0 


(1-w)^ 


I    (v)   e"2iav 


_^ 

v-w 


dv      (n  >  1). 


Since    (E5  ) 


fi       V2      „ 

/(ill)  jjv_   =_rt      (0<w<l), 

J       V  v-w  ' 


iaw  f  <„    T    /-,    ,^V2  1/2 


fo(w)  =  7172— —172  < x  -  rj 

CO 


(l-v)"'     V 

X 


*  ^ 


-iat 

dt  dv 


Employing   (E7) 


/^(l-v)l/2vl/2    .  2  1/2,,       a/2R,         . 

Jo    (v-w)(v-t)        dv  =  -«**  W        (1"W)        6(t"w) 

l/2  l/2 

+rtt        (t-1)  E(t_±)        (0  <  w  <  1,    t   >  0) 

t-w 


26^ 


we  have 


f   (w) 


2e 


iaw 


1/2, .       ,1/2    M 

jtw   '     (1-w  J   '        »- 


■ia    ,    .        l/2^       'i1/2    -iaw 
+  icytw  '     (,1-wJ        e 


+  ia[     ft^Ct-l)1/2         ±\     =  -i«b 


dt 


A  slightly  simpler  form  is  obtained  by  deforming  the  path  of 
integration  into  the  line  joining  1  to  1-i  °°.   Then 


2cce 


ia(w-l)-r-  iti 


(50)   f   (v)   =2ia+- 

jtw        (1-w  J 


5ni    r   (i- 

572"   J  t: 


(^-^)l/2tl/2     e-«*   dt 

w  -    it 


where, as  before,  0  >  arg(l-it)  >  -•=  n. 


With  regard  to  f  substitute  in  (kg)   for  I   from  (^4-1] 


Then 


f>)   = 
n 


oo  . 

-eiaW  ffkZ)1/2   J^      f    ¥JL 

2   1/2, n       0/2    J  k    v  v-w    J         t- 

rt  w  '     (1-w)  '       H)  o 


^  (t)e' 

v 


■2iat 


dt  dv. 


Hence,    using   (E3 ) 


r    i     V2 

+  Jt(tzi)l/2        ^1)      (0<w<1;    t  >0) 

*      4-        '  -f-       _      -i.r  *  ' 


dv  2,l-Wvl/2    _/.        v 

(.v-w  J(,v-t ;  w 


t   -  w 


we   find 


iaw 


fn(W)    =      2   1/2,  1/ 

rt  w  '     (1-w) 


2,l-wxl/2        ,    .      -2iaw 

n    ( )  '      H   (w     e 

v   w  n 


+ 


■27- 


1/2   H  (t)e"2iat 


„  /  %  -law        .  oo      , 

H  (wje  law        p  .  .  1/2 

rtw  '  (1-w)  '    \L 


H  (t)e-2iat 

dt, 


t  -  w 


So  far  the  analysis  is  exact.   Now  consider  what 
happens  when  a  is  large.   The  first  two  terms  of  the  asymptotic 
development  of  the  field  f  are  given  by 


0       ia^w-i— rltl     p     ,-L/^ 

_   /    »        0.  2ae 4        /     t  '    e 

f    (w)^  2ia  +  — tto TTo      /      "1 rz —  dt> 

o 

t  not  being  replaced  by  zero  in  the  denominator  since  1-w  vanishes 
at  the  edge  of  the  disc  w  =  1.   However,  if  0  >  -1  and  v  >  0, 
(52) 

00  00  oo 

p      +0  o-at  p       R     _rif    p      _v(t+i 

'dx  dt 


r  4f?  at .  i  f  t"  e-<*  r  .-*<*«■*>< 

O  O  ^O 

00 

■ivx 


I     7 H5+T       dx 

Jo      (x+a)p 


iav    / 


-ivx 


/  0+1        dx  . 

a         x 

Hence,    taking   0  =  -~,   we  have 


1    . 


i (l-w  )x 


-r-jti  p         -l   l-w 

.   a  eh /  e 

La  +     1/2   1/2  ( n       ,1/2  J  372 

it  '    w  '     (l-w  J  '     va  x  ' 


(53)        f0(w)  —  2ia  +     Wo   i/9         TT7p   / 
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It  is  easily  confirmed  that  this  gives  the  same  distribution 

near  the  edge  as  if  it  were  assumed  that  such  a  region  behaved 

as  a   semi-infinite  plane,  with  edge  along  the  tangent,  placed 

in  the  incident  field.   The  f  and  terms  omitted  in  the  asymptotic 

n 

expression  (53)  represent  corrections  to  this  point  of  view. 

For  f  the  formula  (lj-7)  may  be  employed  in  the  integral  of 
n 

(51).   The  result  is 

Hn(w)e-itXW        h^De-^i*1 

fn(w)~~^ +'  3/2  3/2  1/2 ,.   ,1/2  X 

,    .  .  2  '    it  '  w  '  (1-w)  ' 

(5^J 

-i(l-wjx        h   e  j         -,  /, 

cJ/  rtawJ/  (1-w)  ' 

-3/2 
the  error  being  of  order  ha.   Unless  w  is  near  the  origin 

n 

(J+7)  can  be  employed  in  the  first  term  of  (5k)   with  the  consequence 

(55)        ihn(w)e-iOT       hn(l)e-icW+  b1 

fn(w) 1727; T/2   +  "3/2   3/2  1/2,         T/2  X 

rtw  '  (l-tw)  '       2J'       rtJ/   w  '  (1-w)  ' 

-i(l-wjx        h  .      law 

e  ,       n,o    f   -law       , 

372-  dx  +  — 372  (le    "r7~A/2 

a      xJ/  «awJ/  (1-w)  ' 


The  resemblance  of  the  second  term  of  (55)  to  the  second  term  of  (53) 
will  be  noted.   In  essence,  it  represents  the  effect  of  waves  which 
have  been  multiply  diffracted  (n  times)  at  diametrically  opposite 
points  of  the  edge. 

When  w  is  near  the  origin  it  is  no  longer  to  replace  H (w) 
in  (5*0  by  its  expansion  (k'j)   and  so  (55)  does  not  hold. 
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Suppose  that  ov<  <  1;  then  (kk)   shows  that 

1  °°  i 

v«)-i^2v°»-  ^nl Mil2 

it  x> 

h  _(-ix)e" 


-ax 


^2 2        V0^   dx   +  °(w) 


x      +  w 


1 

=  iwl/2h    (Q)   .    2avei"    f(x£      }l/2     ^ 
it  n^    y  2         J     vl-ix 


P    „3        t /o     h    (-ix) 


2 2       x 

it  ~o  x  +  w 


K   (ax)dx   +  0(w) 


where  l>>ae>>aw.      Because  of  the   small  range   of  integration 

we  can  approximate  K^a  x)  by  l/ax,  hn(-ix)  byhn(0),    e  C     by  land  1-ix  by  1. 

Hence 

/  h*1  r  1/2 

„n(«)K   iAn(o).    Sg      hn(o)/    -^-2    a. 


It  "O  X    +  w 


1  . 


It  X     +    1 

1-  -00  1/2 


•       i  /o  o   V2   4  n  V2 


so  that   H   (w)  -     i  wl/2  h  (0)  {i.  £e  ^        T    x  ^ 

n         ^        n  Jq     x2+1  / 

(56) 

~  _  1  wl/2h    (o) 
~-        it  nv 

xX/2dx         "I 
since   /      —5 =  2     it.      Consequently,    for  a  w  <  <  1, 


r  1/2J 

/  x  '    d) 

'  /      "2 

J  x   +  1 
o 
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(57)  -r*i     P, 

h    (0)       h    (l)e   k     e_2la  h 

„    /    v-u  n  n       n,o 

InWJ- J72      +  03      3/2   1/2   3/2      "      0        1/2 


There  is  no  difficulty  in  confirming  that  the   first 

term  of  fn    when  a  w  <  <  1  cancels  the  last  term  of   f     in    (53) 

1  o 

leaving  a    singularity  at  w   =  0  which  is    smaller  as  a  -» °°  than 

that   in  f    .      Similarly,    the  addition  of  f     produces  a   singularity 

of  smaller  magnitude  than  that   in  f     +  f -,  . 

o    1 

It  is  now  possible  to  obtain  an  explicit  expression  for 
the  distribution  on  the  disc  correct  to  any  order  desired.   All 
that  is  required  is  to  calculate  integrals  asymptotically  to  as 
many  terms  as  are  necessary.   This  calculation  will  be  omitted 
here  since  generally  most  interest  centers  on  various  integrals 
of  the  distribution.   Some  of  these  integrals  will  be  examined 
in  the  following  sections. 

6.   The  High  Frequency  Scattering  Coefficient 

It  is  known  (see,  for  example,  Jones  [29]  )  that  the 

scattering  coefficient  of  the  circular  disc  in  an  incident  plane 

wave  is  -(4/ka  )^A  when  unnormalised  coordinates  are  used  and  A 

is  the  amplitude  of  the  scattered  wave  in  the  direction  of 

propagation  of  the  incident  wave.   Hence,  in  terms  of  our  variables, 

the  scattering  coefficient  a  is  given  by 

1 


(58) 


a     =  -.9  /  wf(w)  dw. 
a.   • ' 


a  Jo 
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Devote  by  o     the  contribution  to   cr  "by  f   .      Then,    from    (5^). 
n  n  '  \s  .>  1  j 


\  ni    p1 


o 


-^[ia+3^2       J      —1/2    J      V2  dXd^ 

it  '  vo      (,1-wJ   '       *"&.     x  ' 


However,    more  terms   in  the  asymptotic   expansion  will  be  required 
and  to  derive  this  we   use    (50)  to  give   further  terms   in  the 
expansion  of  f   .      It   is  a   straight   forward  matter  to   show  that 

-r-Jtl  00  ,  . 

k  p      -i(l-w)x 


It    '      W    '       (l-WJ    '        TL  XJ/ 


(l-w)J_/  Cvtx,        X" 

±  _   3±   +  15  1051  \   dx- 

^        32x2        128x3J 

Wow,    for  x  >  a, 
(59) 

Jo         (l-w)V2  xl/2  ^X        32x2        128x3 

iri 

1/2   ^ 


Hence 


a  Ja,     x  32x^  128xJ 

•(1+s+a-2-J?:3'-  ![r>te] 

32x       J^x3  2X3 


•32- 


-  J    hex  +  a(-  -  -K     +  -At     +    2i_  ) 


W  1+ar  16a1 


-,      .    -ia     -.   -ia 
.   1     ie_         3e -| 

2        2        "        3  J 

a  2a 


-  2.  1  -  ±-  +     J-  cos  a  +     3  +       3   sin  a 

i+a  2aJ  16a  2a  J 


the  error  being  0(a      ). 


The  expression  for     a      (n  >  1)   is,    from   (51), 
(61) 


law 


a      =    £  J    T     f    H   (w)e-laW  dw   +  i-  /U-)1/^1 
n         avLv^nv/  it    J  vl-w 


X 


-2iat 


f    ft-ia/2      Hn(t)e" 


Now 

1 


f    H   (w)e-iaW  dw   =   f  H   (w)e"iaW  dw  -    f  H   (w)e"iaW  dw. 
■J       n  J      n  J       n 


o 


Also 


(62)      /Hn(w)e-iaW  dw  =  g  /Y  (^J1/2  h    (x)e-laX  M(x,w)   dx  dw. 
o  '-„  Jo 

To  evaluate  the   integrals  with  respect   to  w  consider 

V+l„(l)/         v 

z        Hv      (ozj 


2        2 

z  -  y 


dz 


where  -1  <  R(v)  <  ^,   a  >  0,  y  >  0  and  D  is  a  contour  along  the 
entire  real  axis  apart  from  indentations  above  z  =  0,  +  y.   On 


•33- 


account  of  the  asymptotic  behaviour  of  the  Hankel  function  and  the 
range  of  v  the  contour  can  be  deformed  to  infinity  in  the  positive 
imaginary  part  of  the  complex  plane  to  give 

v+l„(l) 


I 


z   H   (      \ 

v  \ajL) 

2   2 

z  -  y 


dz  =  0. 


In  view  of  the  range  of   v  the  contribution  of  the  indentation  about 
the  origin  vanishes  as  its  radius   shrinks  to   zero.      The  contributions 
from  the   simple  poles      z   =  +  y  one  calculated  in  the  usual  way  and 


L 


»  v+1    (1),       , 

x       H        (ax) 

-     dx 


2        2 
x   -    y 


1    .    v/TT(l)/       x  (v+l)rti      (l),        -n±  v 

=  prtiy  (h;    ;(ay)  -   ev  ^(aye*1)} 


c.  rr(l)/      fli  i  -vrti    tj(1)/    \        o   -vrtiT   f    y. 

Since  IT    '(xe      J   =  e  Hx    '(x)   -   2e  J  (xj  we  obtan 


oo      v+1  T    /       > 
r    x       J   (ax)  , 

(63)  /     g   V  g         dx  =  -  |rtyVYv(ay) 

4)  x  -   y 

for   -1  <5l(v)   <  |,    y  >  0,    a  >  0.      In  particular, 


00 


wJ  (aw) 
(6h)  J       -^2__     dw  =  |  «   YQ(ax) 

X   -   w 


.3k- 


and 

law  J 


Pw  JAaw; 
(65)  J    -1T—2  dw  =     2^xYl(aix)' 


O      X    -    w 


The   substitution  of    (6k)  and    (65)   in    (62)   gives,   after  a   use   of  the 
Wronskian  relation, 

00  00 

(66)  ,/Vv).-1**.  -  i/(^)l/2  h^xle-1-  dx  . 

-0  o 

00 

In   /      H   (w)e~laW  dw  the   formula    (U7)   can  be   employed   so  that 

00  00 

I     „   r    \   -iaw    ,  /  /if   x    \l/2        /    %   -iax         n,o        y 

4     Hn(w)e  dv=j  {-(_)/      hn(x)e  -^l/2X 

_(eiax_   le-loocj}    dx   +0(!k)    . 

a 

Combining  this  with    (66)  we  have 

/   Hn(w)e-1-  »   .  i  ^,^hn(x)e--  dx 

+     -S^  /        4/p    (e1CW  -ie  )dx   +  0(-2)  • 

*   a   4        x1/2  a3 


(67) 


The  asymptotic  expansion  of  the  first  term  of  (67)  can  be  deduced 
from  theorems  about  Fourier  transforms  or,  alternatively,  the  contour 
of  integration  may  be  deformed  into  the  negative  imaginary  axis 
and  the  straight  line  joining  1  to  l-i°°.   In  the  latter  case 
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Using  these  results,   together  with   (68),    in   (67)  we  obtain 


rr   t    \  -iaw 
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-^i 


dw 


i  3 
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since  h    can  be  replaced  "by  ih  (0)/8  to  the  order  to  which 
n,o  n   ' 

(69)   is  accurate.   This  completes  the  treatment  of  the  first 
integral  in  (6l). 

With  regard  to  the  second  integral  of  (6l)  we  can  use 
(47)  for  H  (t )  since  only  the  range  t  >  1  is  involved.   Hence 
the  asymptotic  evaluation  of 

1  °°  o  • 

1    p  f  v    J./2      iccw    /k-ia/2   e~2lat  fi  ,   1     a/2,     , ,  * 
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is  required.      Now 
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from  (52).   Consequently  the  terms  in  the  bracket  j   [  of  (70) 
which  do  not  involve  e    J   have  asymptotic  behaviour 
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The  change  of  variable  x  =  a  tan  0  shows  that 
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Hence  our  asymptotic   formula  becomes 
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terms  of  0(h  /a?)   being  neglected. 


The  remaining  term  of  (70)  is 
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and 


1/2 


(t-1) 

t    (t-wj  w 


at 


=  2     {!_!_(!_   w)l/2} 


as  may  be   confirmed  by  the  change   of  variable    z    =    (t-l)        /(l+t         ). 
Hence  our   integral   is 
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terms  of  0(h  /a?)  being   ignored. 


Collecting  together   (69),    (71)  and-   (72)  we  obtain,    from 
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All  that   remains  to  be  done   is  to   find  asymptotic   formulae 
for   h    .      This   can  he  achieved  in  a   comparatively  straight   forward 
manner.      For,   when  n  >  2,    the  combination    of   (43 )   and   (Vf)   gives 
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h 


n-l,o   ^3 
a 


(^  -  1) 


a 

.-^V.d)     9Vl(°)  ,7         . 

13  "  o  l  1" 

32rt2a2  256a  22 


h     i 


a 


-1+2- 


(76)  h  -2ia-rrti  9h        (0)  h 

n  2a  13  n-1  102^a 


2  2  -      ^  2 

8rt   a  a, 


(77)  h     .  -2ia-irti  9h         (0)  h 

n  8a  1|  n-1  20U8q^  g 

8it  a  a, 


The   corresponding  formulae   for  h     are  obtained  from   (32).      Thus 


2      -ia— rrti 


<-     -ia- T~iii  -,-/-i       \  ~o 


al/2(l«) 


so  that 


hx(l)    ^sLe"10^11    +0(a"2), 


2a2 


1 

2        .      1 


.  /-,  >  n        -la-pi     ,    n/      ^  \ 

h'(l)    =  -   — j_  e  k        +  0(a      ), 


W* 


(78)  1  5 

hl(0)    ^e"1-^1    d+|i)    +0(a"    2), 


2 
a 


I  i  -2 

h^O)    =  -   ^  e-ia"   t*1    +0(a"2). 


2 
a 


-43- 


Hence,    from   (73), 


0  .sSsV1^ 


t-jti 


a 


a 


-ia 
e 

2rta 

1 


/   1  .  11     +5i     +Li) 
^    2        32a       16a       8a     ' 


4 


rti 


♦^3<1  +  £  ♦=<§♦§» 


o   2  2 

2rt  a 


•j-rti+ia 

IT 

8rt  a 


(1- 


3i 

4a 


fc  C1  - 1» 


+  0(a        ) 


1  1    . 

92.  "l^1 

+  Li)  +  § 

4a    ;  13 

p^2  2 
2rt  a 


1 

O    O         2  .  1        . 

2$) *      e-ia-^i 
a       1 

2 
a 


-ia  _. 

4na  16a 


■rrti+ia 
h    +i£i)  .  S /         171) 

U   +  I6a;  1  3      U       I6a; 

8rt   a 


_9 
+  0(a  2). 


Similarly, 


2        a 


1    .  1    mj_. 

-■rrti  -rrti+ia 

n22  o   2   2 

2rt  a  oit   a 


h2(0) 


+  0(a  d) 


a 


■ia 


1    . 
2 1  a— rrti 


2rt 


a 


(^l*1!!^) 


-,/-   2      2 
lbrt     a 


■44- 


1  . 


1  5      *l,o 

)    2  2 


-rTti+ia 


h 


i  5         l,o 

n^    2   2 

16*  a 


+  0(a  d)   = 


1    •     K- 

7:    -la-7ltl 

,1  n2e         k 


^ 


X 


1 


X 


•la 


2*a 


^  -2ia— rrti 

{|(1^,+S__!L    j 

2  2 
32*  a 


1    .  1    ... 

j-*i  j-rti+ia 

e  ie 


l  5 

32rt2a2 


?  ? 
128*  a 


+  0(a      )    . 


-!-n 

In   e-eneral   cr     =  0(a  )   for     n  >  1   so  that 

n 


(79)  9 

~2 
a  =  a     +  a     +  a +  0(a     ) 
old 


24) 

a 


,  _,.  .    -ia       _   -1a 

•  «,  1  1  .3i  ,   ie  3e 

ia+  1  -  >-  +  —  2  +  —  3  +  — 5      -3 

4a  l6a  2a  2a 


1  i 

-2ia— r*i  n  •  o2.  .    -ia 


l  3 

1    2  2 

4*  a 


16a  4 


a 


2a 


-2ia— r*i 


2 ,        e  4 


-ia 


/,     ,    49i    ,  1      /n  .   17  iv        e U  "    fi    /-.    p2,        e  4        ,_  S — -  +  — — _ 


1    2  2 

4*  a 


+0(a"2)      =^ 
a 


i      .1      L  i 
ia  +  1  -   77      +  -0   +  — n 
4a     o  2        £   3 
oa       16a 


-2ia— r*i 

e 4_ 

i  3 

1     2   2 

4*   a 


f1   +  Sa  + 


-2ia— r*i     -1 
;  4       \ 


-,,-2  2 
16*  a 


+  0(a  d) 


=  2 


45- 


cos(2a-rit)        7    sin  (2a- -Tit)  , 

1  _  ±      +    -__z —  + ±-     + 


h  2  1  1  11  it* 

m  22  2   2 

Urt  a  6kx  a 


■   j,       sin  (4a  -   ^rtj 

6^ita 


+  0(a"2). 


The   first  three    terms  are   in  agreement  with  the   result   obtained  by 
Levine  and  Wu[l6j;   their  calculations  did  not   extend  beyond  these  terms. 

7.      The   sound-hard  disc 

For  the   sound-hard  disc   on  which  the   normal  derivative   of  the 
field  vanishes,   a   general  expression  for  the  field  is 


u(R)    =  uo(r,z)   +IJ      [u]    ■%-     fiB.,^)  dS 


.er 


where    [u]    =  u(r,    +0)    -   u(r      -   0).      It   has  been  shown   in  an  earlie 
paper    ( Jones  [k])  that   the   integral   equation  resulting  from  the 
imposition  of  the  boundary  condition  oti/cVi  =  0  can  be  converted  into 


1 

cos 


r    ,      2        2>l/2i  1  .    /    i  2  2»l/2"l 

■{a(r     -   x :    ; )   ' '    }  r  d  sinla^  -   x   )       Jdr., 

(ri}  — r-s — 2T172—  toi +  y  g(ri}  d?n  ,  2   2,1" 

(^     -   x    j  ^o  1      (r     -   x    ) 


2    f    (o),    >   v  cosh(a(x  -v    )  '  "} 

~4e  (v)    (x2-/)1/2 


2  ,,2^/2-t 

•dv 


■  k6- 


where   g        (v)    = 


baQ(v,z) 


dz 


and  g(r    )   =    [~u]  .      In  the   earlier  paper 
z=o  1 

we  applied  the  inverse   operator  of  the  first   term.      In  this   case  we 
integrate  the   second  term  by  parts  and  apply  the   edge   condition  g(l)    =  0. 


We  arrive  at 


1      , 


cos 


2 11/2 


(r,) 


,    ,      2        2,1/2. 

ta(ri  -   X    }         \  ■    (n]   sin  a(-/)J 

uTTT^    dri  ■ 1S(0)      (-x2)1/2 

2        2 .1/2 


r    .,      ,      sin{a(r     -   x    )        }  r        (o)  coshfa(x2-y2)l/2 )    . 

,.1/2  1/2  1/2        „ 

Make  the  change  of  variables  r,    =  t  '     ,    x   =  z        ,   v  =  y       .      Then 


r  ,y+i/2 


vL/2 


(-0 


172" 


/(t^)        coslaCt-z)1/2]        dt   _   2.g(o)    sin, 

{     i       ^-^ 

i 

,     f    g(tl/2)        sinfa(t-Z)2)dt  _2         r    (0)(   1/2,      cosh(a(z-t  )^ }   c 

"   XJ      "^^172-  (t_z)l/2  dt    -        zl/2i  r._^l/2 

o 

valid  for  0  <  z   <  1. 


,1/2. 


(z-t)J 


This  is  analogous  to    (6)  but   contains  an    extra    term. 
Nevertheless,   we  attempt  to  follow  the   same  procedure  and  write  the 
equation  as 

r  ktfh  ■-io(t-a)1/L    r  dfeHf)  cosfatt-,)1/2)^ 


W      (t-z)1/2 


dt 


172 


(t-z) 


ITT 


sin  a(-z ) 


•2ig(0)  -  -    l/2 


1/2 


2 
1/2 


:(oV/2)x 


1/2 

cosh{a(z-t )         )      ,, 
(z-t)1^2  ' 


47- 


l/2       l/2 
Multiply  both  sides  by  cos  (a(x-z)    }/(x-z)    and  integrate  with 

respect  to  z  from  0  to  x.   On  account  of  (8),  (9)  and  (B3)  with 

t  =  0  the  equation  becomes 

1/2        1/2  1 

-     Vti/2N        r  -lay  rV+2^ 

g     1/P;        {    /  J-         dy  -  TriH(x-t )  )dt+  /  %±±-r±  rtl   sgn(x-t)dt-2ig(0k 

t  '  -!l/2     1/2     y  4>     t  7 


^L 


t   '    -x 


x2  x  1/2  Z  1/2 

fl^L^dy      =  _2    fcos{a(x-z)        }      r      (o)(tl/2}   cosh  a(z-t)  dtdz> 

_x2        y  4)       z1/2(x-z)1/2   4)  (z-t)1/2 

Instead  of  taking  a  derivative  of  this  equation,   as  we  did  with   (10), 
integrate  the  first  term  by  parts  with  the  consequence,    since  g(l)   =  0,   that 

J        ~7^      {  t1/2    +  XV2  t1/2   -    x1/2      3         =    i  zV2(x.z)l/2X 

o  v 

z  1  /? 

f      (o),   l/2*   cosh  a(z-t)  7         ,,     , 
•    /      g        (t  '     J   s — r-Tp-        dt   dz. 

4)  (z-t)1/2 

2       2 
The  change  of  variables  t  =  w  ,  x  =  v  now  gives 

n  -ia(w+v)    -ia(w-v)       /  ., 

(80)  J    ^   <  ^7-TT-  "  ^TT-)  **  =  g(0)(-) 

o 

where         „ 

(81)  G(o)(v)  =  f      cosja^-z)1/2}   rg(o)(tl/2}  cosh  a(z-t)1/2  ^ 

4)     z    (v  _z)     o  (z-t)  ' 

The  integral  equation  (80 )  bears  a  strong  resemblance  to  (ll)  but 
the  kernels  are  in  fact  different.  However,  an  integral  equation  of  the 
second  kind  can  be  derived  by  similar  means.   Let 

p  -ioc(w-v) 

(82)  G(v)  =  /  g(w)   e  v_v     dw. 


-US- 


Then 

/     v  '     (l-v) 


/    n   -law  1 /     v        (l-v)      "  nf    \    -lotv         ,  C 

S(v)e  =  "      2   1/2,         .1/2  J      ^ G(v)e  dv   +     1/2,.,       ,1/2 

rt    W    '       (1-W  J    '       T3  W    '       (1-W  J 


Choose  C   so  that   g(l)   =  0;   then 

,._,        ,    ,  ,l-w,l/2     eiaW    f  ,   v    ,1/2     G(v)e 

(83)  g(w)  =  -(— )      -g-  /  (i^)      -^7- 


i 

law    p  ,  /0      0/    ,    -iav 


dv. 

'l-v'  v-w 

it      Hd 


An  alternative  choice  of  C  would  be  the  one  which  makes  g 
bounded  at  the  origin.   This  leads  to 


(8.)    *M=-&^   V/<¥>1/2   ^    - 

o 
The  difference  between  the  right  hand  sides  of  (83)  and  (8k)   is 


iaw 

e 

f    G(v)e"iaT 

1  W^Z)^  dv 

2   1/2.         .1/2 
Jt   w  '     (1-w)   ' 

1 

law 

r       1 

""   "      2  1/2,.       ,1/2    v 
it  w  '     (1-w  J 

4    v1/2^)1/2" 

/  ^^^ dx  dv 

/            X    -     V 

■k9- 


from  (82).   This  quantity  vanishes  since 

1 


dv 

"i727  a/ 2,   v"  u- 

r  '  (1-vJ  '  (x-vj 


Hence  the  forumlae  (83)  and  (84)  are  in  agreement  and  either  gives 
g(w)  hounded  and  vanishing  at  the  edge.   Equation  (83)  is  preferred 
"because  it  ensures  that  g(l)  =  0  even  when  approximation  is  used 
and  because  the  behaviour  near  the  edge  is  more  important  than  that 
near  the  centre.   Note  that  (8k)   corresponds  to  the  formula  (lA) 
for  wf(w).   But  it  cannot  be  deduced  that  (83)  would  also  be  valid 
for  wf(w)  by  the  same  reasoning  because  it  is  necessary  that  not 
only  wf(w)  but  also  f(w)  be  bounded  at  the  origin  and  so  (83) 
would  fail  in  this  case. 


With  g(w)  given  by  (83) 


o  00 


G(t)e'iat   _.  , 
— s — *—      dt  dw 

t-w 
■iav  1 


e 

2 
it 


(^2  0(t)e-iat  £;(^)1/2   at 


■50- 


from   (El|).      Hence 

(85,    g(v,  .  -«<°>(v)  -  (i^,1/2  ^7(^>1/2     2i|i£^  at  (ocxi) 


1 

■iav    p  ,  /„        „/.  \   -iat 


which  is  completely  analogous  to  (l8). 

8.   The  Plane  Wave  at  Normal  Incidence 

For  a  plane  wave  at  normal  incidence  g   (v)  =  -ice  and  then, 

from  (81), 

2 

n(o),    ,  0.  r      cos  {q(v2-z)l/2)    ..  .    1/2  . 

Gv   (v)  =  -2i  /    ■-,/g'-'  g- jjo —   sinh  ^    dz 

Jo   z  '  (v  -z)  ' 

=  -   ki   Si(av) 

on  a  use  of  (B3).   Put 

p     -ia(t+v) 
(86)       G(v)  =  -2j     e  t  +y   dt  +  G1(v). 

in  (85).   Then,  because  of  (21), 

^(t+v)  ,fc,/0       -low    P.  „     n/,     -law      _      o   -ia(t-v) 


p      -iat-tv  ,  .      , /0       -iav    p      tt     , /_     -law  n    -iait-wj 

o 
Since    (ElU) 

(87) 


o 
o  o 


r    r  w    f/2  dw  _jt_   ff_v_}l/2_    (_t_)1/2   Hftnn 

JQ    (W  (w"+v)(w-t)      =  t+v      Ul+vj  lt-lj  H^"1^' 

00  .1  „    /    \   -iaw 

„   ,    .        9,Ut.1/2     -iav/%  t    ,1/2  er10*,,        ,1-^1/2  e'^ffW    }l/2  Vw)e         dv 
Gx(v)   =  2(~)  e         J    (-j)  ^     dt   -    (— )  -        j    (^)  w+v  dw 


-51- 
Putting 

i4v  a/ 


)  G   (v)    =   (i±£)^   e"laVM(v) 

1  v 

we   obtain  -, 

(89)       M(v)   =   i±(v)   -  -   /      (^)  —  dw 


where 


VT>  ■  2.   fe»      t» 


,  /0       -iat 


dt 

1 


(90) 

1  00  i 

=  ^7%^    dt 

J        t1/2(l-tv  -   it) 


by  deforming  the  contour  into  the  straight  line  joining  1  to  l-i°°. 

Equation  (89)  is  analogous  to  (3O).   To  find  a  solution  valid 
at  high  frequencies  we  try 

00  ,  V 

u,f    \        1   v  ^/2   f  Ll      -2ia(t-v)  ,.  ,  M  /  » 
M(v)  =  (^)     J   ^7  e  dt  +M2(v). 

o 

Then 

v  ,,„  p  L(t)e-21a(t-T)    T  1/2 

o 


T    /,  v   -2icct  «>  T    /,  n   -2iat 

Li  (t)e  p      .      ,/„  L   (t  Je 


T5T"       «-./%' 


J        t+v  J      t-1  t+v 

o  1 

1  /    ,  -2: 

1    f  ,1+  vq/2     _2 


M0(w)e 


(iTJl)^     _£ dw 

It  J         1-     W  W     +   V 


o 


-52- 


where  ad-vantage  has  "been  taken  of   (87).      Choose  L     so  that 

00  »        * 

P  -2±ct{t-v )  -2iat  ,         1 /0 

r  L  (t)  { 1 f^-)dt  =  (— )V2f(v). 

./  1  *■    t    -    V  t    +  V   J  V  1 


This  integral  equation  is   studied  in  Appendix  D   (DI7)  and  its   solution 
(D2U)   is 

00 

Lx(y)   =  g^Y(±T)l/2      ^(x)e-iaX  M2(x,y)  dx 


o 


where 


M  (2)(ax)Jn(ay)   -   4x2H   (2)(ax)   J  (ay) 
(91)     MQ(x,y)   = 


^    '*'   ~  2      2 

x  -y 


Hence  in  this  case  we  form  the     iteration    scheme 

00  t    (+  \   -2ia(t-v) 
,  ;       p    L   (,t  Je 

(92)  M   (v)   =   (T^-r/2   /     -^ at, 

w  n  1+v  /      t   -   v 

(93)  Ln(y)   =  ^f^^)1/2   in(x)e-lcK  M2(x,y)   dx      (n  >  l), 

°°  L        ft)   "2iat 

fel0   ,n(v),.|(-iI,V2^i_e      dt  („>2) 

and  i,    is  given  by   (90). 

It  may  be   shown  in  a   similar  way  to  that   of  Section  h  that 

1  °°  -\  /o 

O 

(95) 


x  IC(ox)J  (cty)  -   xyK   (ax)J  (ay)-i 
{— -± ^-p p-S —  '    dx. 


2  2 

x     +  y 


■53- 


From  this   it   can  be  deduced  that,    if  ay  is  not   small, 

(96)     Ln(y)-    ¥if}^  «bW   +    "  1  . 

(32rtaJ  y 


9.      The   Distribution  on  the   Sound- hard  Disc 


Let  the  contribution  of  G(v)-G   (v)  to  g  be  g     and  the 

contribution  of  M    be  g   .      Then,    from   (83)  and   (86), 
n  n 

1               /  °°          /         \ 

,        ,  /0     law      p            1/2  -iav    p    -ia(t-v) 

(97)  «  Or)  =  2(^)1/2  £,      /  (JL)  ^—    /  f—          dt  dv 

o                     w                   2      J      jj7  v - wj    t   -   v 


o  o 


and  from  (88) 


roftl  t    )  A-^l/2  eia      f    ,Hv ,1/2       Mn(v)e-2iaVH 

(98)  gn(v)   =  -(— )  — ^      (— )  -^  dv. 

Wow,    f or  0  <  w  <  1  and  t   >  0,    (E12) 

1 

/nn,       f    ,  v    xl/2  dv  2/  w    vl/2  .  /.        »  /   t    a/2   H(t-l) 

(99)  /     U — )         / n — tv  =  rt  (1 — )        6  (t-w)  +  Tt[-—T)  '     — -. 

J         1-v  (v-w)(v-t)  1-w  t-1  w   -   t 

o 

Therefore 


Also 


■5k- 


1  °°  T    /,  s   -2iot 

law     p  ,  /0     ,  p    L   (,t  Je 


o 

oo  . 

(101)  =  Ln(w)e  -    (— )  e         J  (^^)  * 

n  1 


L    (t)  0.     , 

n  -2iat,. 
e  dt. 

t    -  w 


We   shall  not  attempt  to  take  these   forumlae   further  although 
approximations   for  large  a  on  the  lines  of  Section  5  are  possible. 
Instead  we   shall  concentrate  on  the   scattering  coefficient. 


10.    The   Scattering  Coefficient   of  the   Sound-hard  Disc 

1 


The   scattering  coefficient   of  the   sound-hard  disc   is   -2^  /g(w)  wdw. 

Denote  by  cr     the  contribution  of  g   .      Then,    from   (100), 
n  °n  '  ; 


■ictt 
a    =  -  sm  I 


« 


r   o.         -     1/2H       ^/2     iaw    f  ,   t    ,l/2  e'1"^^    . 

{-2w   +  -  w  '     (1-w)  '       e  /    (— —)  ' dt}  dw 

fl  J      t-1  t-w 

1-   \/2   /     w^d-w)1/2   e4        /  ^  (i.i+^)   dx  dw 

it  '     4)  J      x  '  ^2x 


x  '  32x 

a 


Now,    for  x  >  a, 


■55- 


1/2,-,       a/2      iwx 


f    w1/2(l_w) 


e         dw 


■      3„.  13. 

1/2  1X-Urtl  ,.         1K  2^ 

(1   +  ^  +—  2}    +  ~172    (1 
4X        32x  2xJ/ 


2x 


3/2 


2i) 


Hence 
(102) 


a   ~  2# 
o 


.      1 
32x  HX        32x  X 


2K 


-IX 


X  OX  X 


2ft 


-,        .  -,  .    -ia  -ia 

,     ,    . /I    ,   1         ,      1    1        ie  e 

a  a 


4a  oa 


1  1        .  ^  1 

1  -   — 5     -   ~2   sin  a  +  -_      cos   a 

^a  a  a 


Turning  now  to   cr   (n  >  l)  we  have 


(IO3 )  1 


0^,2^ 


{wL    (w)e  dw  -  "(l-w) 


,      iow    n  /      \l/2  L    (t)   -2iot 
1/2ri_T7a/2   e  /    /_t_\  j^_  e 

it         J    \t-l)  t-w 


dt}dw 


Wow 

(iou)   n1 


T     /    v   -iow  t    /    \   -iaw  /  /    ,    -law    , 

wl^wje  dw   =    /   wL    (,wje  dw  -    /   wL   (,wje  dw. 

1  o 


Also 


■56- 


00  00 


wLn(w)   e"iaW  dw  =  ^  /  w   /    (i^)1/2in(x)e-iax  |^(x,v)  dx  dw. 


00  /  . 

^    w  J   (aw J 


— ^ 2     dw  =  ^jtYfax)  and    / 

x      -   w  4) 


w   J   (aw; 

2        2 

x  -   w 


dw 


I^Cax) 


the  right   handside  becomes,   after  a  use   of  the  Wronskian  relation, 


00 

1      f'      1/2 ,.       .1/2  „     /    v   -lax   . 
-  — .    /      x  '     (1+xJ  '      I     (.xje  dx 

111 4) 


e   h     I    (0) 
n 

"    0   1/2.      3/2    ' 
2rt  '     i   a 

In  the   first   integral  of    (101+ )  the   form   (96)   can  he  used  and,    furthermore, 
the  asymptotic   expression   (2/itoy)  cos(ay  -  r^)   for   J,  (ay)   can  he 

employed.   Hence 

1    . 

00  00  ,      .     J"3*1 

TwL    (w)e-iaWdw~    [(iw^d^J^i     (w)e"iaW   +  -^%     cos(aw-^)    }   dw 
J        n  ^      rt  n  ^awi/2 


22i      -23Ti"ia 
Jta 


£  (i)  -  -    (h  (D  +  £«(i)} 

n  a       4  n  n 


1   . 

V°)e  3 

-^ sin  (a  -   -^t). 


kn 


a 


Consequently 
(105) 

n  ■  o1/2  -ice 

r    r    \  -iaw,      ~  2       e 

-    /   wL   (w)e  dw  ~  

,/         n  «a 


^-lwi)  +  £,J^ 


a  (o)  r1 

g sm(a-jj;rtj 

1+rt 


L 


a 
1    . 

— rJtl 

e   *     in(0) 
2n1/2a3/2 
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In  the  second  integral  of  (l03 )  the  form  (96)  can  be  used  with 
J  replaced  by  the  first  term  of  its  asymptotic  expansion.   We  thus 
obtain 

(106) 

1  °° 

1    f    1/2 /_       sl/2   iaw    P ,   t    a/2   e"   latr   l,t+l\l/2.    , .  v 
*        W        (1"W)        6         ./    fe}  t^  v    {   ^"H        'n(t)    " 


I      (0)  O-      4- 

^—,h    (1   +ie2lat)}   dt   dw. 
ckat^ 


Consider   firstly  the  terms   inside  the  brace    (      }  which  do  not   involve 

2ictt 
e  .      Wow 


1     i  1 

2/.,       x2      . 


1  ,11 

3  °°  —  /  J_0O—  — 

^/n       xc      .  fiti  p      2/.     .     a/2   -ax  titip      2,  .2      ia-ax 

w   (1-w)        iawn  4      /     x   (1-ix)  '    e  ,  4      /     x   (1+ix)     e 

— * —         e       dw  =  e        /      — !l— - — dx  -   e        /     —\ — : — L dx 

t-w  J  t   -   ix  a        t-l-ix 


TTti+iat     p      -itx  i  _    . 

u       x 

a 

1  °°        /         \ 

^i+iat    p     i(t-l)x     (  , 

J^       (5!+S-2-)dx 

a 

from   (52).      Hence  the  terms   in  the    {      }   of   (106)  not   involving  e        '  have 
asymptotic  behaviour 


4x  ^   (t-l)1/2     *  8rtcct 


_  1/2/  3/2' 

2ir  '      J  a     x  ' 


-la    p      -ix  .  _. 

•  e  /  e  n  ,.    ix       31  /,    .    ix,  •> 

~2~J  "~3/2    {1  ie 

a 


Ux* 


d-^^x)-^)^ 


'•it  n 


B^  ]  ,  .a/2 

(x+aj 


21/2  ! 

2(x+a)J/    it      4     n  n 


dx 
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The   substitution  x   =  a  tan      Q  shows  that 


dx 


cc 


(x-ttx)   '    xJ/ 


(2-21/2), 


a 


a 


(x^)1/2x5/2       =        2 


dx 


2      ,3 


^)17r7^   =  ?   2V2  ■    j" 


Hence  the  value  of  the  integral  is 


(107) 


e 

rta 


■la 


l/p  ^n^1^    7  -  ii5n(l)     -  1/? 

<2l/2- 2>  V1' +  -t-  <5  -J/a' +  -S-  f|  - 2 ' ) 


-ia 


a 


The  remaining  term  of    (106)   is 
(108)      j^    f  „V2(l.w)l/2eiaw| 

1 


8rt  a       o 


dt 


t(t-l)1/2(t-w) 


dw. 


The   inner  integral  is  a   function  of  w  with  no   singularities   in    (0,l)  ]_ 


except  near  w  =  1  where  it   equals    (l-w) 


.       ..       ,-1/2     f      dt 
Is    (l-w)      '      J      -57; 


tJ"/"(t+l) 


+  0(1)   =  rt(l-w)    +  0(1). 


Hence  the  first  term  in  the  asymptotic     expansion  of   (108)   is 
(109) 


in(0)ela 


8rt 


a 
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Combining    (105)^(107)  and    (IO9)  we  have 


(HO)  .  .  -o1/2,   Ml         -r«i 

r--ia  -7.  -ia  _.  i2   '     jj    (0)  h 


CT„    ~     2ft 


£_    {(2-Ii-|_    )  j  (1)  -  |i  i-(i)  + — jr-^—}  -  e    l/o  ,lo  i  (0 

rta      l        oa       2a  nx  2a     nN  8a  n   1/2   3/2      nv 

2it       a 


0/2 

The  error  in  this  expression  is   0(i  /a        ). 

It   is  also  convenient  to  have  available  asymptotic   formulae   for   2,   , 

n 

From   (9^)  and    (96),    for  n  >  2, 

0    (    \  f  I   t    V1/2      e"2l0Ctfi    fl+t,l/2  /    »  Vl(°j      ,  2iatu., 

vv)  ~  -  j  fei}     t*  {* (— }     Vi(t)  -  ~^;t3/2 (1+ie     )}dt- 


It  is  sufficient  for  our  purposes  to  observe  that  this  implies  that 

1  1 

-2io6-riti     22rt2^     .  (l)  ii     .  (0) 


/    n  1  "4  n-1  n-1  /  dt 

(2a)1/2(l+v)    '       8lta       '    {        t(t-l)1/2(tfv) 


e 

n  '    '  it 


o  /p  p 

for  v  >  0,   the   error  being  0(i        /a        ).      The   substitution  t    =  1    +  x     shows 

that 

00  00 

/     T7o =  *(l  -  --I  /o)>   that    /  ~p T7o   =  o  n  and  that 

{     t(t-l)1/2(t+l)  21/2  J±  t2(t-l)1/2 


/  T7p p  =  ^ s/?)-      Hence 

-{  t(t-l)1/2(t+l)^  25/2 

^ti"2ia  U     .(0) 

^     "^  -Turv-2    Vi(1)+^    ^^1/2^ 

2it  '    a  2  ' 

?,1-2Ia  ii     ,(0) 

(112)       ,__  (1,    „    e-^      in_i(l)  .  1-i_  (1      5     ,, 


i+TC 


a 
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(113)     i  (0) 

n 


■j-rti-2ia 


1/2  1/2 
it  '    a  ' 


n-1 


ii     -,(0) 


The   corresponding  formulae   for   f,    are  deduced  directly  from   (9O) 


and 


11.. 
2  "T^1"^ 


a    ^     ^  ^n   e 


2*"e  i(l-v)    , 


so  that 


*]_(!)    = 


1 

2 
it   e 


P   -j-rti-ia 


a 


1/2 


{1-K)(ct      )}, 


1/2  -T^i-ia 

t        e  ,.,  .^/   -1. 


ii(l)    =  "   ^172  (l^(cc"X)}, 


2a 

1/2   -tJti"ia 
.    /  ~  \        ail       e 

*1(0)   =  —I72 


[1   +0(a"1)). 


a 


Hence,    from   (110), 

11.. 
o   — r  Jti-ia 

a 


1  1    . 

r--ia  .  -ia       n2.  "It  1 

£         r2_  i  •  _  £-       +  Li] 
Jt  ^     8a       2a  V 


1/2   1/2 
it       a 


0(a"3), 


a2   =  2/& 


-ia 


Tpci 


^       "2(1)   "   L^2a3/2      "2C0) 


+  0(a" 


2e 


ia       J*1"21"^!)       1^(0)  1 


=  2*  ^     {~  Z^Krr^ 


2n    '    a 


-210,8.(1) 

"(l"21/2))+L1/2a3/2jtl/2al/2 


+0(a"3) 


1     1 


2  -jfiti-ia 


a 


fiti-3ia 


.    -ia  ,                     -2ia~ 

e ,   le /,  1        s 

~172        +  ~1J2    (1  "   "1/2 > 

itJ/            2na  '  2  '             2nd 


1/2 


+  0(a"3) 


from   (ill)  and   (113).      Similarly 
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a     =  2Kf-       l   (1)   +0(a-3) 
3  ita  3 

T-rti-2ia 

=  ^"7172  1/2       -1        i2(D+0(a-3) 

2ir       a 


1   •    c> 

pri-bia 


■2Kh7^575-  +  0(a"3)- 

2it        a 


1     , 
-p-n-1 

In  general     cr     =  0(<x  )   so  that 


(1U) 


a  = 


a     +  a 
o 


1   +  a2   +  a     +  0(a"3) 


=  2R 


-  ^rri-Uia  1    .    „.  -prti-2ia 

,     .  i        1        +  ie2  21  ^     IT1-21-  3e^ 

"      1/2   3/2   S  +  o   1/2  5/2 

n       a  on  '    a 


a 


i+a         «a 


ie 


rrti-6ia 


2*3/2a5/2 


+0(a"3) 


=  2 


1 1i2  0/2   sin(2a-^it)   +  -g   sin(ikx-   g-ir) 2 

rt  '    a  rta  *kx 


+  aJ2cp/2°3{Za  '  *°  ^^T2?^  sln^~h 


+  0(a"3) 


This  formula  is  in  agreement  with  the  result  obtained  by  Levine  and 

Wu[l6]  except  that  their  factor  of  the  term  a     cos (2a  -  rt)  is  7 

3 
instead  of  75.   The  cause  of  the  discrepancy  has  not  been  traced  but 

is  very  likely  due  to  an  arithmetical  error.   Dr.  Levine  does  not  think 

it  possible  that  the  next  stage  of  the  approximation  adopted  by  Levine 

and  Wu  would  produce  a  term  of  comparable  order. 
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11.   Uniqueness 

In  this  section  we  show  that  the  solution  of  the  integral 
equations  (l8)  is  unique.   This  immediately  implies  that  the  solution 
of  (3O)  is  unique.   Furthermore,  it  has  been  already  shown  that 
(ll)  implies  (l8)  i.e.  a  solution  of  (ll)  provides  a  solution  of  (l8). 
On  the  other  hand,  if  F  is  any  bounded  solution  of  (l8)  it  is  easy 
to  verify  that  f  defined  by  (lk)   satisfies  (ll).   Therefore,  there 
is  an  equivalence  between  the  two  integral  equation  for  the  types  of 
solution  under  consideration.   Hence  a  unique  solution  of  (l8)  implies 
a  unique  solution  of  (ll).   But  since  (l)  possesses  a  unique  solution 
(Jones  [27]  )  which  is  a  solution  of  (ll)  it  follows  that  our  process 
provides  the  unique  solution  sought. 

We  turn  now  to  the  proof  that  the  solution  of  (l8)  is  unique. 

Since  f(w)  is  continuous  in  0  <  w  <  1,  (13)  shows  that  F(v)  is 

bounded  and  continuous  in  0  <  v  <  1.   Near  v  =  0 

e 

F(v)  =  f(0)  f   SSU0(1) 

o 
where  1  >  >  e  >  >  v.   Hence 

F(v)  =  f(0)v  in  (e/v)  +  0(l) 

so  that  F  is  bounded  at  v  =  0.   Near  w  =  1,  f(w)  =  A(l-w)"2  +  0((l-w)    ) 

and  so  1 

1       _± 

F(v)  =  A^    %^}  2dw  +  0(1)  =  0(1) 

1   §  1_e 

since  /   - — 775- =  n,   which  differs  from  the  given  integral  by 

(l-w)1/2(w-v) 
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0(l).   Consequently  a  solution  of  (l8)  is  required  in  which  F(v)  is 
continuous  in  0  <  v  <  1. 

Suppose  there  are  two  such  solutions  of  (l8).   Then  their 
difference,  F  (v),  is  of  the  same  type  and  satisfies 


(115) 


FQ(v) 


2  -iav 
v  e 


1      \ 
(1-w) 


tf(l+v)1/2  4)  w1/2(w+v) 


Fo(w)e_iaW  dw. 


for  v  in    (0,l).      It   follows  that 


F   (v) 
o 


1 

■1 


V 


-      it(l+v) 


1/2  mxlFol 


(i-wr 
i/2/     ; 

w  '     (w+v  J 


dw 


<    {1  - 


1/2 


V 


If  max   IF   (v)|    occurs  at  v  =  v     we  have 
1    o        '  o 

Y    1/2 

mi|Pj<    tx--a_yi/sj 


max  F 


If  max    |F   |    =  0  there   is  nothing  to  prove  and  if  max    |F    |    ^  0  we  must   have 


v     =0.      Now  as  v  -»  0 

1 

2     p      F   (w)dw 

Fo<'»  =  -J,/     -TI2T        +o(l) 


o        w        (w+v) 
e 


1/2 

-        F    (.) 

j:         o   s 


dw 


"b  w  '     (w+v) 


+  o(l) 
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where  0  <  £  <  e,   by  the  mean  value  theorem.   Since 
'-b  w  '  (w+v)    v  '  v 


(116)   Fq(v)  +  Fq(£)  =  o(l) 

from  which  may  be  deduced  that   F   (0)   =  0.      Therefore  max  |F    I  =  0  and 
u  o  '    o ' 

F      =      0  i.e.    the   solution  of    (l8)   is   unique, 
o 

A   similar  method  may  be  employed  to  determine  the  error  caused 

by  stopping  our   iteration  process  at    some  n.    Let   F   (v)   =  v  '     e  I    (v)/(l+v) 

Then  ,  -, 

—  1  - 

v  h    (v)e                   2   -iav              r    (1-w)   F   (wje 
^   /    \               n                         v  e                      /                     n 
Vv)  =  " TTo  '  — ~T73        /      TTT. : dw- 


o 


(l+v)1/2  "    rt(l+v)1/2       4)  w1/2(w+v) 


As  above 


1  1  11 

|F   (v)|    <  2  2|h    (v)|    v2    +   {1  -   v2/(l+v)    }  max|F    |. 

If    I F   |    has   its  maximum  at  v  =  v     /  0, 
1    n1  n 

_i  i  l 

max|F   |    <  2  2|h   (v)|  (l+v    f  =  0(a  "n+  2) 
1    n     -  '    n    A  '  n 

from   (76),    (74)  and   (78).  On  the  other  hand  if   |  F   |    has  its  maximum  at  v   =  0, 

then,  as  v  -»  0,  F  (v)   +  F  (£)   =  o(l)  by  an  argument   similar  to  that   leading  to 
n  n 

(ll6)  and  F     =   0.      In  either  case  we   can  say  that   the  error  in   stopping  the 

_n+l 
calculation  at   I     ,    is  0(cc~        2)   in  F. 
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The  integral  equation  (85)  requires  a  more  elaborate  treatment 

and  only  the  proof  of  uniqueness  for  \oc\    <  it  will  be  given  here.   Firstly 

the  possibility  that  G  has  a  logarithmic  singularity  at  the  origin  cannot 

be  excluded.   The  difference  between  two  possible  solutions  G  satisfies 

o 


tj.     n  /o       ~±av     r       4.       1  In       G    (t)e"1 

(117)  <.„(,)  .  -(i?)1/2  I  /(^     _o_  dt) 


which  may  be  written  as 


(118)    Go(T)  .  .  (i?)i  f  ""/(i)1  ^(tje-^fi  -  ^,  at. 


The  form  of  the  solution  needed  is  G  (v)  =  A  In  v  +  /(v)  where  ?(v) 
is  bounded  and  continuous  in  0  <  v  <  1.   To  examine  whether  such  a  solution 
of  (ll8)  is  possible  determine  the  behaviour  of  the  right  hand  side  of 
(118)  as  v  -»  0.   Wow 

p    vG  (t)  -iot  p      vG  (t)e"lat 

4  /ihL?i%*)  dt  =4  Ai-^'2^)  dt  +  0(v)- 


In  the  integral  over  (0, e)  use  the  assumed  form  of  G  .   Note  that 
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V7„(t)e  7„(lJe  P 


dt 


P    v7Q(,tJe  7QlUe  p 

I  t^^i-t)1/2^) dt  '  a-e)1/2      v  -1    J^WJ 


%«>*'**     l/2         .x      r 

=  TTo v  tan         V~ 

(1-|)1/2 


v 


=  0(v2). 


Further 


1  1/1 

p6        n"2  n^r£/V     n"2      ,      , 

/    vt      mt  dt  =  v        /       t — ic£vtl  dt 

4)      t  +  v  4>         t+i 


ni     e/v       -i  -^  -i  n  "I 

=  vn  2J       ft1""2-   tn"'2+,,,+(-)n-1t"2+  ^^      )    ^n(vt)dt 


1 

■-• 


2(-J     v  /  ^ dt    +  0(vj 

H)  t    +1 


=    (-)n  Jt  v       '   ^n  v   +  0(vn+1   ^n  v)    +  0(v) 
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00 


since    /  —x —       dt    =  0.      Hence 

J  t    +  1 

o 

(119)  1  _       /,     X  .        ,  1 

o        e  dt  „    .      1/2  ,   2> 

-it  Av'       in  v  +  0(v    ) 


■-      t 
o 


i7W[7Vv) 


l/2  2 

ind  the   0(v  '     )  term  is  actually  rt7o(0)v   +  0(vin  v).      Consequently 

rG(t)e-iat 
GQ(v)   =AMv-T2      /        l/2     ■       l/2     ^    +  0(1). 

JTV    '  -D       t     '       (1-t  )    ' 


One  concludes  that  G  (v)  can  satisfy  the  given  conditions  only  if 


o 


(120)   r1  G  ft)*"1"* 

J  ~M^F2  dt  =  ° 

o 


and  it  then  follows  from  (ll8)  that  G  also  satisfies 

o 


1 


o1  G   (t)r-i0Ct 


2  §  o    G   (tje 

(121)  G   (v)   =  J      (l+v)2   e-1<XV   /        °        l/g   l/g dt. 

-0      (l-t)1/2t1/2(t-Hr) 


Since    (l2l)  can  be  rewritten  as 

Go(v)   =  _(1+v)   e  /    (__)        ___ {  __     +  __}   dt 


it  follows  from  (120)  that  G  must  also  satisfy 


1    .     1     1  r    /.  s   - 
,2     e"laV  T/l-t.2   Go(t)e 


iat 


(122)  o0w  -  (^r  5       <^r  -^-     dt. 
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The  facts  that  G  must  comply  with  the  condition  (l20)  and 
o 

also  be  a  solution  of  (122)  are  consistent  with  the  discussion  in 

§  7  of  "the  relation  between  (83)  and  (8U).   But  there  a  certain  form 

for  G  was  employed  whereas  here  (120)  is  a  consequence  of  (117)  and 

the  general  conditions  imposed  on  G. 

It  will  be  noted  that  (122)  differs  from  (115 )  only  in  the 

sign  of  the  right  hand  side,  but  this  difference  is  sufficient  to 

prevent  the  application  of  the  method,  used  for  (115),  to  (122). 

It  is  still  true  that  max  I G  I  occurs  at  v  =  0  because  if  G  contains 

1  o  '  o 

a  logarithmic  term  it  is  infinite  there  and  if  it  is  bounded  the  same 
method  as  for  (115)  may  be  employed.   However  (ll6)  takes  the  form 
7(v)  -  ?(£,)  =    o(l)  and  we  cannot  deduce  that  max  |G  |  =  0. 

If  we  multiply  (117)  by  G*(v)e2laV  v(l-v2  j"1'2,  the  star 
indicating  a  complex  conjugate,  and  integrate  with  respect  to  v  from 
0  to  1  we  obtain 


1  11  1        ,    ,     *,    ■>   -iot+iov 

r  v 1- ,  ..2  2iav^       1  rr(     vt     ,2  Go(t)GQ  (v)e 


rt  J 
o 


^lZt)l/2G0(t)e"iafc-xtdt!2dx 


Since  the  right  hand  side  is  real 

1 
(123)  /     — 2  x/2      |Gq(v)|2   sin  2av  dv   =  0. 

If  |a|  <  t«  this  is  possible  only  if  G  (v)  =0  so  that  uniqueness  is 
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proved  for  |a|  <  pit. 

An  extension  of  this  range  can  be  obtained  by  observing  that 

«   GQ(v)eiaV  p1     GQ(t)e-iatdt 

^  -,  /0   -,  /„      /       1  /0   -,  y0 is   integrable  at  the  origin  on 

(1-v)1/ V/2    4)        (1-t)1/    t1/2(t+v) 

1 

account   of   (119).      Hence  if   (l2l)  be  multiplied  by  G  *(v)   e2l<XV/v(l-v2)2 
and  the   imaginary  part  be  integrated  from  0  to  1  we  obtain 


r  lG  (v)l2  •    o 

/,„,  %       /      '    o     _J sin  2  ay 

lt)  (,1-v     J 

A  combination  of  (123)  and  (124)  is 


^  lGn^^!2        2      .2 

°  -  Wo   (i4.it  v  -  — )  sin  2avdv  =  0 
(l-v2)1/2  V 


which  implies  G  (v)  =  0  for  |a|<  it. 
o  '  '  - 
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APPENDIX  A 


In  this  appendix  we  derive  the  formula  (3)  for 
n  exp[-ia(r  +r   -  2rr  cos  (cp-cp-.)}  '      J 

/        — YF dcpr 

J  [r  +r      -2rr     cos    (cp-cp   ) }    ' 

Suppose  r  >  r     and  write  "b   =  r  /r   so  that  b  <  1.      Then  one  integral  we 
have  to  calculate   is 

pU     cosfyfl+b   -2b   cos    (cp-cp   ) }    ' 

/         " 5 173 dcPi 

J  r[l+b   -2b  cos    (cp-cp   )}' 

where  7   =  CCr.      Put  cp   -cp  =  9  and  then  in  the   interval    (-cp,0)    replace   9 
by  2n+9.      The  result   is 

2jt  „  ..  I 

P        cos!r(l+b    -2b   cos   9)    '    ] 

J 5 TB —    de 

0  r(l+b^-2b   cos   9)    ' 

p2*     cos(7(l+b2-2b   cos   2X)1/2} 

=     k  /  - dX 

J  r(l+b    -2b   cos  2X) 

by  changing  9  to  2jt-9  in  the   interval   (n,2it)    and  then  putting  9=2X. 

The   substituion     ~ —  =  tan  9  now  gives 

cos   2X  -  b 

P      cos[7f(l-b2sin29)1/2-b   cos   9}] 

~  J        2 2~l72 dG 

0  r(l-b   sin~9)    ' 

o2       cos    {7(l-b   sin  9)   '    }cos(7  b   cos  9) 
=  k  /  g 2~T72 dQ 

£l  cos    (a(r2-x2)l/2}cos(a(r2-x2)l/2} 

=  TJ  ,    2      2,1/2,    2      2,1/2  dX- 

^0  (r   -x   )    /     (rx-x   )    ' 


-  71  - 


When  r  <  r   just  interchange  r  and  r   in  this  last  formula  because  of 
the  invar iance  of  the  original  integral  under  such  an  interchange.   Hence 

p*  cos[o;{r  +r  -2rr  cos  (cp-cp  )  ]  ' 

/ T~2 172 —  dcpi 

J  [r   +r  -2rr  cos  (cp-cp  )}  ' 

min(r>rl)  Cos(a(r2-x2)l/2}cos{a(r2-x2)l/2} 

=  kJ  f    2   2,1/2,  2  2,1/2  CbC- 

^0  (r  -x  )  '    (rx-x  ) 

The  integral  involving  the  sine  can  "be  dealt  with  similarly 
but  for  our  purposes  it  is  more  convenient  to  proceed  otherwise.   We 


have 


r-      sin[a[r  +r   -2rr  cos    (cp-cp   )}    '  J 

/        2~~2 172 dCpl 

J  (r  +r   -2rr   cos    (cp-cp,)}    ' 

/2     P2rt  J^frtr^-^cos    (cP-cp1)}IL/2J 

=    (o  ««)    '       \  '    g     g ± JTJJ dcpx 

J  {r  +r   -2rr     cos    (cp-cp   )}    ' 


2n       |  it 


=  aJ   J 

0        0 


?      ?  l/2  ~~ 1 

sin  OJ    \a (r  +r   -2rr     cos    (cp-cp   )}    '      sin  OJdOdcp-!^ 


after  using  Sonine's  first  finite  integral  for  Bessel  functions    (Watson 
[30]).      Use  of  the   addition  theorem  for  Bessel  functions  now  gives 


1 
,2* 


2«a    /  sin  0  JQ(a  r   sin  0)J0(a  r1  sin  Q)d9. 

0 

Now  it  has  been  proved  elsewhere  (Jones  [k] )  that 
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1 

[Z   sin  9)  sinV  U  0  cos   GdO 


2 
/     J    (z  sin  9)  J    (I 


v-—  u+- 

2  ^2 


v-|   X+|   J*  Jm((z2^2coB28)^2) 

=  z   d  Z  d   /    J    (Z  sin  9)  -£± 

0  2  2   2  t^ 

(z  -Z  cos  0) 


.  u-X+1  .     2X  _,,, 
sm      9  cos   9d9 


where  $(X)  >  -  '-,   $0)  ><#(x)  -  1,4V)  >  -  |.   Using  this  result  with 
v  =  —,      |i  =  -  — ,  X  =  0,  z  =  Qlr  and  Z  =  ar,  we  obtain 

-*         ./„        .      ^    .       ,..,   2     2  2  Ml/? 


>2       cos(ar   sin  9)sin    [a(r   -r     cos"  9)    '    ) 

71 — 2     2q>iA ^ de- 

^  (r      -  r  cos   9)    ' 


Hence 


Qtt  O       Q  "I   /o 

p      sin[a{r  +r   -2rr  cos    (cp-cp   ) }    '""J 

7      ,  2+  2  0 ; — -TT2 — dcPi 

wn    (r  +r  -2rr  cos  (cp-cp,  ) }  ' 

r       ,2  2vl/2   .    ,  2  2,1/2 
n      cos  a(r  -x  )  '   sm  a(r  -x  )  ' 

=  kJ  ,   2  2x1/2,  2  2x1/2  dX- 

0    (r  -x  )  '  (r^-x  )  ' 

Since  r  and  r  can  be  interchanged  on  the  left  hand  side  without 

affecting  the  result  the  same  must  be  true  of  the  right  hand  side.   This 

can  be  confirmed  in  the  above  procedure  by  taking  z  =  cc  r}      Z  =  en  r 

An  independent  verification  is  to  note  that 

{f   2  2xl/2  ,  2  .2x1/2. 
n        exp  ((z  -a  )  '  -(z  -b  )  7  } 

J       /  2  2x1/2,  2  .2x1/2      dZ  =  °- 
°_     (z  -a  )  '    (z  -b  )  ' 
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Here  a  and  b  are  real  with  0  <  a  <  b  and  the  branch  lines  are  (-«>, -b), 

2  2  l/2 
(-00, -a), (a, 00)  and  (b,<»)  with  (z  -a  )  '   =  ia  at  z  =  0.   The  contour  of 

integration  lies  just  above  the  real  axis  and  the  equation  follows  by 

deforming  the  contour  to  infinity.   If  we  take  the  imaginary  part  of 

this  equation  we  find 

P    .    .   2  2.1/2     ,.2   2xl/2 

/   sm(a  -x  )  '      cos  (b  -x  )  ' 

J  1    2  2x1/2,.  2  2x1/2  " 

£     (a  -x  )  '    (b  -x  )  ' 

n  ,  2  2,1/2   .   ,_2  2xl/2 

/   cos  (a  -x  )  '   sm  (b  -x  )  ' 

{    (a^x2)1/2  (b^x2)1/2 
and  the  restriction  a  <  b  can  obviously  be  dropped  now. 
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APPENDIX  B 


In  this  appendix  we  consider  the  evaluation  of 

r         t„t       ^/2,  -ia(t-z)1/2 
/  cos_fa(x-z)  '    \e 

/  —   i/4 —  1/5 dz 

J0    (x-z)1/2(t-z)1/2 

where  x  and  t  are  non-negative  real  but,  for  the  purposes  of  this  cal- 
culation, a  can  be  an  arbitrary  complex  number.   Write  the  integral  as 

x  l/2        1/2  1/2        l/9 

r   exp(io:(x-z)  ;  -ia(t-z)   }  +  exp{-ia(x-z)  '  -ia(t-z)  '  } 

1/2      1/2 
In  the  first  integrand  make  the  change  of  variable  (t-z)    -(x-z)    =  y 

1/2      1/2 
and  in  the  second  integrand  put  (t-z)   +(x-z)    =  y.   Thus  we  obtain 

n(t-x)1'2   -lay       ntt-x)1/2  -iay 

7  -r- dy  -  J  ~Ay- 

tl/2.xl/2  t1/2^1/2 

If  t  >  x  the  integration  in  both  is  along  the  real  axis  and  does  not 
include  the  origin.   If  t  <  x  the  upper  limit  of  integration  is  positive 
imaginary  (with  the  convention  on  square  roots  adopted  in  the  main  text) 
and  the  lower  limits  are  negative  and  positive  real  respectively.   In  this 
case  the  contour  of  integration  can  be  deformed  into  the  real  axis  with 
an  indentation  round  the  simple  pole  at  the  origin.   Hence 

x  ,  1/2        tl/2+xl/2 

F   cosfafx  z)l/2)e-ia(t-z)  F  e'iay 

(Bl)    /  cosfa(x-z)    )e   dz  =   /         e_^  dy  _  niH(x_t) 

JQ  (x-z)1/2(t-z)1/2  V/2_xl/2    y 

the  principal  value  of  the  integral  being  understood  and  H(x)  being  the 
Heaviside  unit  function. 
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Since  (Bl)  is  valid  for  arbitrary  complex  a.  we  can  replace 
a  by  -a  and  add  the  resulting  equation  to  (Bl).   In  this  way  we  obtain 

tl/2+  1/2 

(B2)  f   cos(a(x-z)l/2}cos(a(t-z)l/2}  dz  =   T        cos_av.  ^.^h^ 
J0         U-z)1/2(t-z)1/2  t^-x1/2    Y 

Similarly 

,1/2^  1/2 

X  1/2  l/2 

(B3)  f   cosfa(x-z)   )sin(a(t-z)   }  dz  =  f         sin  ay  dy 

J0      (x-z)l/2(t-z)l/2  Jt1/2-xl/2    y 

Interchanging  x  and  t  in  (B2)  gives 

.  l/2_,_.l/2 

1/2  1/2  X         * 

(B4)    T     co*W*-*)  /g)°osfa(x-z)       )   dz  =      r  co^a^  H(t_x)- 

J0  (t-z)1/2(x-z)1/2  JxV2_tl/2  y 

The   subtraction  of    (B^)    from   (B2)    leads  to 

n      cos(a(x-z)1/2}cos{Q;(t-z)1'2} 

(B5)       YJ2 172 dz  =  ~ni  sgn  (x_t)- 

(x-z)  '  (t-z)  ' 


t 
since 


x^-t1/2 
f  £OS_av.  dy  =  Q> 

t^-x1/2 


-  76  - 


APPENDIX  C 

This  appendix  is  concerned  with  the  solution  of  (ll)  at  low 
frequencies  when  a  is  small.   The  method  of  approximation  is  that  of 
truncation  of  the  series  expansion  for  the  kernel.   The  error  so  intro- 
duced can  be  estimated  but  will  not  be  considered  further  here.   To 

2 
illustrate  the  process  we  neglect  the  powers  above  a     in  the  kernel. 

The  resulting  integral  equation  is 

f     wf(w)f-J^  -  2ia  -  a^\    dw  =  p(°)(v) 
^       Iw  -v  > 

or         x 

(CI)      f   Z^1   dw  =  F<°>  (v)  +  a0a  +  a/ 

0  W  "V  ! 

where  the  constants  a  and  a  are  defined  by  aQ  =  2i  /  wf(w)dw, 

1  0 

f      2  2  2 

a     =    /     w  f(w)dw.        Making  the   change   of  variables  w     =  x,   v     =  y 

0 
in   (Cl)   we  have 


fiM„.^ 


2 
+  a.joc  +  a.,  a 


x-y  0  1 

0 

which  is  of  the  same  type  as  (13)  with  solution 

xi/2f(xi/2)  = i f  zl/2(i-y)l/2  fF(o)(yi/2} 


{F(0)(yl/2)  +  ^   +  a^2}  ^ 


-li- 


on  choosing  c  so  as  to  eliminate  the  singularity  at  x  =   0.   Since 

r    (1-v)1/2 


$?- 


o  y  (y-x) 


2 
a^a+ano; 


'   (I-,)1/2 


i/2   "er^*      i     f   (i-y)x/%(Q)f  i/2, , 

f (X    }  -   ,_   .1/2  -   2,n   ,1/2  J    T7T^ 
jr(l-x)  '     it   (1-x)  '   °_   y  ' 


Thus 


2 
a^a+a^a: 


cr  "i"  2  r    d-v")'^  -,(Q) 


*    fi-v2)1/2 


(C2)   f(v)  =  —I rr/2"  -   2n   2,1/2  J    2  2  F   »dv- 

jt(1-w  )  '     it   (1-w  )    n   v  -w 

2 
The  constants  a„  and  a  can  be  determined  by  multiplying  by  w  and  w 

and  integrating  over  (0,1)  in  the  two  cases. 

For  instance,  with  an  incident  plane  wave,  we  have  from  (19 ) 

F^  (v)  =  k   cos  a  v  +  h  a  v  Si(av). 
In  principle,  there  is  no  necessity  to  approximate  F    in  (C2)  but, 
for  simplicity,  we  shall  use  only  the  first  two  terms  in  its  expansion  in 
powers  of  a   i.e. 

F(0)(v)^  k   +  2a2 v2. 


Since 


and 


J  2     2    dV  =  "  2  * 

0   v  "w 


f    (1-V2)1/2  V2  ,      1   ,1     2, 

J     22 dv  =  -«(--  v  ; 

0    v  "W 


we  have 


-  78  - 


2 
aQa  +  a^  k  1  2      ^ 


(C3)      f(v)  =  ;   2.i/2  +   .    2,1/2  t1  "  I  °  (I  "  W  )] 
jt(l-w  J  '     jt(l-w  )  ' 


Hence 


1 


21  /  wf  (w)dw  =  —  (a^a  +  ana2)  +  —  (1  +  \-  a2) 

,  i  Jt      U       X         Jt  X^i 


ao    J 

0 
and 

1 

P        2         1  2        12 

a  =  /   w  f  (w)dw  =  r  (aQa  +  a  a  )  +  1  +  g  a  . 

0 

The  solution  of  these  simultaneous  equations  is 

ao  =  f  [1  +  ^  } 

On  it 

a1  =  l 
correct  to  the  terms   in  a  and  a     respectively.      Substitution  in   (C3)    gives 

f(v)  Jfl  +  ^  +  (i.M  o?ni-v2r1/2  - ^r  ^1/2  +  °^3)- 

Jt  Jt  \d  d.i  Jt 

Jt 

This  is  in  agreement  with  the  results  of  Bazer  and  Brown  [7]  . 

2 
The  calculation  of  f  correct  to  0(a  )  involves  the  determination 

of  an  and  a  the  first  two  moments  of  f .   To  obtain  the  terms  which  are 

0(a)  it  would  be  necessary  to  solve  three  equations  for  the  first  three 

moments.   Clearly,  although  all  the  integrals  required  are  relatively 

simple,  this  approach  does  not  provide  the  expansion  in  powers  of  a   as 

rapidly  as  the  integral  equation  found  earlier.  (Jones  \_kj  ) . 

Finally  note  that  the  determination  of  a„  is  virtually  equivalent 

to  finding  the  scattering  coefficient. 
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APPENDIX  D 
In  this  Appendix  the  solution  of  the  integral  equation 

n  f     -ia(w+v)    -ia(w-v)Ni 

(Dl)      /  f(w)  < + >dw  =  h(v)       (v  >  0) 

K  J  I   w+v         w-v    J 

0 

is  derived.   We  find  the  solution  such  that  f  is  integrable.   The  method 

of  solution  is  not  difficult  once  certain  integrals  containing  Bessel 

functions  are  available.   These  we  shall  now  obtain. 

Let  ix  p   /p\ 

Uxy  H^;  (ax)  J  (a  y)  -  hY    v\   >    (ax)  j  (ay) 

(D2)      M(x,y)  =  ■ 2~~2 

x  -y 

where  the  usual  notation  of  Bessel  and  Hankel  functions  is  employed. 

Consider 

n  ia(z+w) 

r  /  M(z,y) dz  where  y  >  0,  w  >  0  and  C  is  the 

k    J  z+w 

C 
contour  shown  in  Figure  1.    The  only  singularity  inside  the  contour 

is  the  simple  pole  z  =  y.   Hence 

p       gia(z+w)       2  gia(y+w) 
k   J    M(z'y)  T^ dz  = 


a   y+w 
C 

after  using  the  Wronskian  relation 

HJ2)(ay)  JQ(ay)  -  H^2)(ay)  J^ay)  =  2i/itay. 

On  the  large  circle  of  radius  r  the  asymptotic  formula 
(2)      (z   iV2  -i(«z-  |  v«-  |n) 
Hv  ^■)/°V«j    e 

(valid  for  -2rt  <  arg  z  <  it)  may  be  employed  and  then  we  see  that  the 
contribution  of  the  large  circle  tends  to  zero  as  r  ^».   On  the  straight 
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Figure  1 
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portions  of  the  negative  real  axis  we  put  z  =  x  e—   and  use  the  fact 

ufhxe711)    -  Eff  }(xe-rtl)  =  -2i  sin  vn   Hf  }(x)  +  lie™1  Jy(x). 
For  the  indentations  about  z  =  -y  and  z  =  -w  we  need 
H(2)(xe'Ti)  +  H(2)(Xe"rtl)  =  2  cos  v*  H(2)(x). 

There  is  no  contribution  from  the  indentation  around  the  origin  as  its 
radius  shrinks  to  zero.   The  net  result  is 

P        -ia(x-w)       r     ia(y+w)    -ia(y-w)  ^ 
(D3)      7  K(x,y)  — dx  =  §(  — +  S_. j  _  |  niM(w,y) 

0  y 

where 

hxy  J  (ax)  Jn(ay)  -  %  Jn(ax) j  (ay) 
(DlO      K(x,y)  = ^— 


2  2 

x  -y 


Similarly,  it  can  be  shown  that 


n  -ia(x+w)        r    ia(y-w)    -ia(y+wh 

(D5)      /  K(x,y)  £— dx  =  %   { +  2__ \  -   ±  jtiM(w,y) 

J     '     x+w         a    I  y-w        y+w     J   2    v  ' 

0  J 

■o  n       •        ti(2)         ,      iccz  ,       TI(l)  -iCte  ,,     . 

By  replacing  H         and  e  by  H         and  e  we   obtain  m  a 

v  v 

similar  manner 


o                    ia(x-w)  r  -ia(y+v)        ia(y-w)^ 

(D6)  j    K(x,y)   S_ ^  .  tl^—  +   ^— }♦  |  ^(v.y), 


P  ia(x+w)  r  -ia(y-w)         ia(y+w)  ^| 

(D7)  /    K  (x.y) ; dx  =  -    < +  — >  +  ±  JtiK.  (w,y) 

v    |;  J         v    'J;        x+w  a    ^     y-w  y+w  I        2  lv    '"" 
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where  M  is  the  same  as  M  with  IP   taking  the  place  of  IP 
Adding  (D3),  (D5),  (D6)  and  (D7)  we  find 


(D8)    /         <    cos  q(x+w)+  cos  a(x-v)  \  ^  =  h  j   cos  a(y+w)  + 
0 


K(x,y)  {  ^^  *+  ■  ^3^—  >te       a  y+v 


cos  Qi(y-w) 
y-v 


«L(w,y) 


where 


(D9)      L(x,y) 


hxy  Y1(0Qi)J0(ay)-ky   YQ(ax)  J-^oy) 


2  2 
x  -y 


so  that 

(D10)     M(x,y)  =  K(x,y)  -  iL(x,y) 

and 

M1(x,y)  =  K(x,y)  +  iL(x,y). 

Next, it  follows  from  (D8)  that 

-ia(w+v)     -iQl(w-v) 


(DID 


i  r  fe-ia^+v> 
dv  =  -    /  2—z + 

w+v  w-v      J        n     J      \        w+v 

0  V- 


L(v,y)  I*- 


-ia(w-v)' 


w-v 


4  f  cos  a (y+v)   cos  a(y-v) 
a  \   y+v         y-v 


K(x,y) 


cos  a(x+v)   cos  a(x-v) 
x+v         x-v 


dx 


dv. 


!,  W 


"  -ia(w+v)    -ia(w-v) 
e  e 


w+v 


w-v 


cos  a(y+v)   cos  cc(y-v) 
y+v         y-v 


dv 
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1  f     fe-ia^+v)  ,  e-^lfcosa^)  +  cos  g(y-v)\  dy 

2  ,/   1   w+v        w-v    f  \   y+v         y-v     / 


and  the  integral  on  the  right  hand  side  can  be  evaluated  by  expressing  the 
cosines  as  exponentials  and  using  formulae  such  as 

00 

F        e       V  dv  2R,        v      2iay   ,       .  ,?iav  2iayw,        , 

/        (Y-w)  (v-y)    =   *  °(v-^    e  +   rtl<e  "  e  }/(w-y)- 


Thus 


(D12) 


-    fe-ia(v+v)  +  e-ia(w-v)>     f  cos  a(y+v)  +  cosa(H|  ^ 
/      1       w+v  w-v         J    )         y+v  y-v         J 

2R,        ,             f sin  a  (w+y)        sin  «(y-w) 
=   Jt  5(w-y)    +   rt    {   :7r± —        "  777T 


^        w+y  y-w 

for  the  values   of  y  and  w  involved.        Inserting  this  result   in   (Dll)   we 
obtain 

p  f  -jcc(w+v)  -ia(w-v)  ~^  |  ,  (    .         /    .    \ 

/  _ ,        v    /   e      v  ,    e     _J I    ,  4it      ,        n    ,  i|/  sin  a  (w+y) 

L(v,y)        —— +  — - — >  dv  =  —   &(w-y)   +  ±{  


w+v  w-v 


■    r   av  =  —    hiw-yj    -t-  — \   

J  a  OK     J '      a\^      w+y 


00 


(m3)  _   Bina(y-w)  |      _rtK(Vjy)    _J    K(x?y)    I  sin+«(x+w)    _   sin_a(x-w)  1  ^ 

=  -^  b(w-y)    -itK(w,y) 
from   (D3),    (D5),    (D6)    and   (D7).      Since, from  (D5)  and   (D6)  , 

P                       f     -ia(w+v)             -ia(w-v)  "^ 
J     K(^)        [~-^ +  ^ )dv  =   -iK(w,y) 

we   see,   recalling   (D10) ,   that 
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p  C   -ia(w+v)    -ia(w-v)>\        , 

(mk)  /  M(v,y)    *—— —   +  ^— — -  V  dv  =  -  ^  o(w-y)  . 


w+v       w-v  a 

~0 


Multiply  (Dl)  by  iaM(v,y)/l*jt  and  integrate  with  respect  to  v 
from  0  to  oo.  on  account  of  (D1*0 

00  oo 

/  f(w)5(w-y)dv  =  g   /  h(v)M(v,y)dv. 
0  0 

Hence  the  solution  of  the  integral  equation  is 

oo 

(D15)     f(y)  =  jg  J    h(v)M(v,y)dv. 

0 

p  p 
When  en  ->  0,   cnM(x,y)  ->  8iy/:rr(x  -y  )  and  in  the  limit  (D15)  gives 

CO 

(B16)     f(y)  --%     \      ^f- 

n      o   v  "y 

That  this  is  indeed  the  solution  with  integrable  f  and  a  =  0  can  be  confirmed 
by  an  independent  analysis.   The  integral  equation  is  now 


/   §£Ma»-h(v). 

0 
As  in  Appendix  C  the  general  solution  of  this  is  given  by 


W  V2x       1     f  vl/2h(vl/2)  ,   +   C 
f  (v   )  =  -  -£-572  j   ^-^ dv  +  -^ 


or 


f(w)  = 


Jt  W        n  w 


f  4m  du  +  c_  _ 


2    J      2   2      w 

JT  W    °„      U  -W 


On  choosing  C  so  as  to  eliminate  the  non- integrable  singularity  at  the 
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origin  we  recover  (Dl6).   The  presence  of  the  constant  C  in  the  general 
solution  here  indicates  that  in  suitable  circumstances  there  may  be  so- 
lutions of  (Dl)  in  which  f  is  not  integrable  but  the  left  hand  side  has 
a  meaning  in  spite  of  that. 

For  the  related  integral  equation 


00 


r   -ia(w-v)    -ia(w+v)  A 
(D17)     J     tWlS— £__ |dw  =  h(v)      (v>0) 

0       ^ 
start  from 

(D18)     M  (x,y)  =  ^ 1— £ ± 2 . 

x  -y 

By  the  same  method  as  the  preceding  we  show  that 

P  -ia(x-w)        f       ia(y+w)    -ia(y-w)N) 

(D19)      /   K  (x,y) dx  =  -{-S— - S  _  ireiM  (w,y), 

J    lv  ,J/   x-w         a  (^   y+w         y-w   J       2   2V  ,J  ' 

0 


"       -ia(x+w)       r  ia(y-w)   -ia(y+w)  ^ 
(D20)    J  ^(x.y)  2-_ dx  =  |  {§__  .  e__ I  +  _  ^^  ? 

0  V 

o         ia(x-w)       r  -ia(y+w)    ia(y-w)^ 

(D21)     /  K  (x,y) dx  =  ^  p—r +  %   jtiM_(w,y), 

v        J   lv  '      x-w       a  I   y+w        y-w    I   2    3 

0  J 

P  ia(x+w)        r  -ia(y-w)    ia(y+w) -\ 

(D22)     J  Ki(x,y)  ^__  dx  =  g  |  ^—  -  ^-  ]  -  g  rti^y) 

0 

where 

UxyJ  (ax)J  (ay)  -  ^x2j  (ax)J  (ay) 

(D23)     ^(x^y)  =  ^ L__ 1 2 

x  -y 
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and  M  is  the  same  as  M  with  H    replaced  by  IT   .   It  follows  from 
(D19)-(D22)  that 


(D24)    f    K   (x,y)  {  cos  a(x+v)  -  cos  a(x"v)  I  cbc  =  M  COS  a^< 
v       J        1K    'J '  \         x+w         x-w     J  a\        y-w 


0 


cob  a(y^?\  +    ( 

y+w      J  lv  ' 


7) 


where  K^x^y)  -  iL1(x,y)  =  M2(x,y). 
Since 


00 

e 


-ia(w+v)    -ia(w-v)  \     C  ,        y  ,        .~\ 

v      e K    J  cos  «(y+v)   cos  a(y-v)  \. 

w+v       w-v     |  1    y+v        y-v     J 


0 


=  A(w-y)  -  rtf sln  a^  +  slna(w+y)l 

1^       w-y  w+y         J 

for  w  >  0,   y  >  0  it  follows  from    (D2*0    that 

p  f*   -ia(w+v)  -icufw-v)^  ,  ,,  (  .        .       s 

f     L,(T,y)|s i 1    .£__ lldv  =  ^  S(w_y)_^/^^__(Z_Z1  + 

J         lv    '•'M         w+v  w-v  I  a     v     ■"     a\     w-y 


sin  a(w+y)  1  rr  ,        v         f  T,  /        \    J    sin  a(w-x)    ,    sina(w+x)l  , 

+    w4  ,y/) +  rtVv^  -j  Vx>y)  ( -w — +    w+i  7 to 

0 

=  "^  5(w-y)  +  itKj^Cw.y) 
from  (D19)-(D22).   Also 

p  f   -ia(w+v)    -ia(w-v)  ^ 

J  V^r  w+v     - e  w-v    j  dv  =  «"i^ 

0 
so  that 
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f        r        i  f  e-ia(v+v)    e-lQ:(v-v)\  .      U«l  _,    . 
J  M2(v,y)  ^    v+v    "    v.v   j  ^  =  "  —  B(v-y). 

0 
Hence,  if  (D17)  be  multiplied  by  M  (v,y)  and  integrated  from  0  to  <x>, 
we  obtain 


f(w)  ^j-  &(w-y)dw  =  /  h(v)M2(v,y)dv. 
0  0 

Consequently  the  solution  of  (D17)  is 

00 

(D2U)     f(y)  =  ]£-   f  h(v)  M2(v,y)dv. 

0 

It  is  of  interest  to  note  that  in  the  course  of  the  preceding 

analysis  it  has  been  demonstrated  that  the  equation 

P             r    -ia(w+v)    -ia(w-v)  ^ 
f (v)  =  X  \     f(w)  I    — + >■  dw 

0       V 

has  an  eigenvalue  X  =  i/n  and  a  corresponding  eigenfunction  f (v)  =  K(v,y) 
The  integral  equation 


-ia(w-v)   ^-ia(w+v) 


f (v)  =  X  f      f(w)  {  -  V  dw 

0 


has  an  eigenvalue  X  =  i/jt  with  associated  eigenfunction  f  (v)  =  K,  (v,y) 
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APPENDIX  E 


In  this  Appendix  an  alternative  evaluation  of 

J        \--J  (v-w)(v-t)  (0<w<l) 


0 
is  given.  A  Fourier  transform  of  the  integral  is 

■let    r    /i-v\l/2 

0 


iPt    /    f^l      T^tT  « 


/   [  i=2      "i  e      sgn  P  dv 
\  v  /        v-w 


Take  branch  lines  from  -oo  to  0  and  from  1  to  «>.   Then  we  specify  the 

l/2         1/2 
branches  of  v    and  (l-v)  '   by  saying  that  |  arg  v|  <  jr  and  that 

(1-v)  '2   =  (v-1)1'2  e^-1'2'*1  with  2it  >  arg(v-l)  >  0.   Wow,  for  P  >  0, 

extend  the  interval  of  integration  to  (-00,00),  the  path  lying  entirely 

below  the  branch  lines.   Then  the  right  hand  side  is 

0O  ,  y  OO  /  00  / 

r    /i      \1/2      .      -iPv  n,      ,\l/2  -i0v  o,-,,     ,1/2         ipv 

r  (in  dv  +  A—     £LJL—  av  -  A—    -^h^—  av. 

J   \  V  /       V-W  J  \  V  /      v-w         J  \  v  /      v+w 

-co  10 

When  the  path  of  the  first  integral  is  completed  by  an  indentation  below 
v  =  w  the  contour  can  be  deformed  to  infinity  to  give  no  contribution 
so  that  the  value  of  the  first  integral  comes  from  the  indentation. 
Consequently  we  obtain  when  (3  >  0 

1  0 
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When  p  <  0  we  extend  the  interval  to  (-00  j0o)  but  with  the  path 
lying  above  the  branch  lines.  On  repeating  the  process  and  remembering 
the  presence  of  the  factor  sgn  p  we  find  (El)  also  when  P  <  0. 

An  inverse  Fourier  transform  now  supplies  the  original  integral  by 

Pfl-v]1/2  dv        1_   f   ipt  f   2  fWl1/2   -i(3v 

J  \v  J  (v-w)(v-t)  =  2*  J      S    {*     U)     6 

0  -00      v 


^a/^-ipv  ?,       .1/2    ip,     >> 

I— J    v~w- dv  -  *V  w     v^r dv)  de- 

1  0 


Since     —    /     e1^     dp  =  &(t)   we   obtain,   for  0  <  w  <  1, 
—00 

0 

,1/2 


H(tzll+      /2>t\   ^Sj^tl 

t-w  V  -t  /  t-w 


In  particular, 


to)         /  (^)1/2  T^TT^T  =   4^)l/25(t-w)   ♦   „(^)1/2  Still 


v  /  (v-w)(v-t)  \  w  /  \   t  /  t-w 

0 

(0  <  w  <  1,   t  >  0), 


^m    r  fi^)1/2    ay  /i±t^1/2 

l^       J       \vj         (v-w)(v+t)      "  t+w     \  t  / 


(0  <  w  <  1,   t  >  0), 

\      v     y  ^  v  — w  y  ^  v  r  u  /  u~rw         \     u     / 


By  multiplying   (E^t-)  by  t  and  letting  t  ->»  we  find 
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m  f  (i^)l/2  &j  .   -„  ■  (0  <„<!). 


Then,   for  0  <  w  <  1, 

}     ,,    a/2  1/2  r1/-,,    a/2    .  ..     ,1/2        .        ^ 

<*)      /     (£)U        "-/{^)        ^+t(^)      T^M" 

0  0     ^ 


2   l/2h       >1/2R/+      u      ,1/2/,     nq/2   H(t-l) 
+  k  v  '    (1-w)   '    &(t-w)+jt  t  '    (t-1)   '      — ; L- 


t-w 

from  (E2)  and  (E5).   In  particular 


1   (l-v)1/^1/2  „       +   2  1/2,,   ,1/2,.   .  +  ^(t-l)1/*    ^ 
•S — ^-T7 tt—  dv  =  -ir  +  it  w  '  (1-w)  '  5 (t-w)  +  — r * ^   H(t-1, 


(E7)   J  7;iw)(v-t)  dv  =  -rt  +  *  * '  d-^  '^c*-*) +    tV    ;    hc*-1) 

0 

(0  <  w  <  1,  t  >  0), 

p  ,_  ,1/2  1/2         ,1/2,,  _q/2 

0 
By  interchanging w  and  t  in  (E2)  we  obtain  a  formula  for  the 
integral  on  the  left  hand  side  valid  for  all  w  and  0  <  t  <  1.   Thus 
the  only  range  not  covered  is  that  in  which  neither  w  nor  t  lies  in 
(0,1).   Elementary  computations  show  that  for  these  w  and  t 

r  f^)1/2     av      -*-  f(^\1/2  t^)1/2\ 

J     {   v  J  (v-w)(v-t)  ~  w-t  H  w  J        "  V  t  J        J  . 

0 

Inspection  of  this  and  (E2)  reveals  this  can  be  conveniently  combined 

with  (E2)  to  give,  if  w  /  (0,1) 
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(E9: 


m 


1/2 


dv 


(v-w)(v-t)  "  w-t 


w-l^2   A-lW2 


w 


(H(t-l)+H(-t)} 


Thus,  for  w  >  0  and  t  >  0 


(E10) 


dv 


7    (v+w)  (v-t)  '"  t+w  \l  t 


¥)  Wfe) 


1/2 


It  may  also  be  shown,  in  a  similar  way  to  that  used  for  (E6)  ,that 


(Ell) 


(1-v) 


1/2  l/2 


(v+w) (v-t' 


dv  = 


0 


^  [w1/2(i+w)1/2+t1/2(t-i)1/2H(t-i; 


for  w  >  0  and  t  >  0. 

By  changing  the  variable  v  to  1-v  we  deduce  that 


(E12) 


v  l  '       dv 


1-v/    (v-w)(v-t) 


n 


i^f2^-^  +  ^k{£if2mt-i)n(-t: 


(0  <  w  <  i) 


(E13) 


t-w 


.  \V2  (    t  ,1/2 


w-1/     Vt-1 


H(t-l)+H(-t)} 


(w  /  (0,i; 


and,  in  particular,  that 


(Elk) 


1/2 


dv     =  _n_ 
1-v/     (v+w) (v-t)  "  t+w  \   \l+v)         '  \t-l 


1/ 


'  t  W2 


H(t-l) 


(w  >  0,  t  >  0) 
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